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Promoted by the recent progress of Berry phase physics in spin galvanomagnetic communities, 
we develop a systematic derivation of the reduced Keldysh equation (RKE) which captures the 
low-energy dynamics of quasi-particles constrained within doubly degenerate bands forming a single 
Fermi surface. The derivation begins with the Keldysh equation for a quite general multiple-band 
interacting Fermi systems, which is originally an Ni^ by Nt, matrix-formed integral (or infinite-order 
differential) equation, with Ni, being the total number of bands. To derive the RKE for quasi-particle 
on a Fermi surface in question, we project out the fully occupied/empty band degrees of freedom 
perturbatively in the gradient expansion, whose coupling constant measures how a system is disequi- 
librated. As for the electron-electron interactions, however, we only employ the so-called adiabatic 
assumption of the Fermi liquid theory, so that the electron correlations effect onto the adiabatic 
transport of quasi-particles, i.e. the hermitian (real) part of the self-energy, is taken into account 
in an unbiased manner. The RKE thus derived becomes an SU{2) covariant difi^erential equation 
and treats the spin and charge degrees of freedom on an equal footing. Namely, the quasi-particle 
spin precessions due to the non-abelian gauge fields are automatically encoded into its covariant 
derivatives. When further solved in favor of spectral functions, this covariant differential equation 
also suggests that quasi-particles on a doubly degenerate Fermi surface acquire spin-selective Berry 
curvature corrections under the applied electromagnetic fields. This theoretical observation gives 
us some hints of possible experimental methodology for measuring the SU{2) Berry's curvatures 
by spin-resolved photoemission experiments. Due to the non-trivial frequency dependence of (the 
hermitian part of) self energy, our RKE is composed of Berry's curvatures in the d-\-l dual space, 
i.e. k-oj space, so that the dual electric field is already introduced. To provide a simple way to 
understand this "temporal" component of the (7(1) Berry's curvature, we also provide the dual ana- 
logue of the Ampere law, where the "spatial" rotation of the electric field in combination with the 
"temporal" derivative of the well-known magnetic component are determined by the 17(1) magnetic 
monopole "current". 

PACS numbers: 71.10.Ay, 71.27.-|-a, 79.60.-i, 72.15. Gd 



I. INTRODUCTION 

Gauge fields often appear in effective low energy the- 
ories in condensed matter physics, whenever the sys- 
tem's low energy manifold is restricted by some "lo- 
cal" constraints^. Classic examples include doped Mott- 
Hubbard insulators, two dimensional electron gases in the 
fractional quantum Hall regime, low-dimensional quan- 
tum spin systems and highly frustrated magnets. 

In these "strong-coupling" problems, the constraints 
are implemented locally in real space. For example, the 
on-site Coulomb interaction in the Hubbard model or 
t-J model forbids double occupancy on every site. In 
strong magnetic fields, the wavcfunction for electrons in 
a two dimensional electron gas must be annihilated by 
the Landau level lowering operator, which is a differen- 
tial operator defined locally for eacii value of the two- 
dimensional coordinate. In the quantum dimer model, 
spins are presumed to be bound into singlets which cover 
each site exactly once. In many cases, these local con- 
straints take the form of a version of "Gauss' law" im- 
posed on suitably introduced "electromagnetic field" and 
"charge" variables. Due to the Gauss' law constraint, 
the low energy effective theory for such systems becomes 
mathematically equivalent to a form of (usually compact) 
quantum electro-dynamicsii^ii^. 



Another kind of constrained system occurs in the weak 
coupling region, such as in a Fermi liquid (FL), where 
the states at the fermi surface are supposed to be well 
separated from the other fully occupied/empty bands by 
a sufficiently large (direct) band gap. The low energy 
manifold is thus spanned only by the Bloch states of those 
conduction bands forming this Fermi surface. Even in 
the presence of strong electron correlations, we may still 
assume that the low energy Hilbert space is spanned only 
by those quasi-particle excitations constrained within a 
Fermi surface, as far as electron-electron interactions can 
be introduced adiabatically in comparison with the direct 
band gap. 

In most of the literature, Fermi liquids are assumed 
to be well described, by ignoring electron correlations 
completely. In such a non-interacting case, the projec- 
tion to a single low-energy band or degenerate bands has 
been studied extensiveljs^iii^i^ii^, and yields an effective 
equation of motion (EOM) for the conduction electrons 
in this band, say the a-th band, moving under the in- 
fluence of external electric fleld e and magnetic field b. 
Such an EOM contains an dual magnetic field , which 
acts on a quasi-particle like a Lorentz force in k-space: 
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z is a CP''^""'^ vector in FLs with A^-fold degenerate con- 
duction bands, i.e. z = (zi, Z2, • • ■ , zat) with z^z = 1. 
This complex-valued vector describes the internal de- 
grees of freedom associated with the degeneracy at each 
fc-point. Correspondingly, the dual magnetic field, mag- 
netic gauge field A"' and "orbital magnetization" Mq, 
are all N hy N Hermitian matrices. In a non-interacting 
fermi gas, they are defined solely in terms of the peri- 
odic part of Bloch wavefunctions of conduction bands 
\ua-) (cr = 1, • • • , iV £ a); 
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with Ea being the energy dispersion for the a-th band. 
Especially, the third EOM dictates that the CP^~^ vec- 
tor z precesses due to the nontrivial matrix structures 
of [Ma] (Zeeman field) and [A"] (Wilczek-Zee phase) re- 
spectively. 

One purpose of this paper is to enlarge the regime of 
validity of this effective EOM, so as to include metals for 
which electron-electron interactions are significant, i.e. 
Fermi liquids rather than Fermi gases. The SU{N) ef- 
fective EOM mentioned above is clearly valid only for 
a non-interacting Fermi gas, or within a mean-field de- 
scription of ordered states such as (ferro)magnetic metals 
in which the ordered moment does not fluctuate at all. 
The fundamental framework of Fermi liquid theory, how- 
ever, implies that these effective EOMs should be prop- 
erly generalized into a realistic metal, where the electron- 
electron interaction/magnetic fluctuations are not weak. 
In fact, Haldane has recently argued that the "renormal- 
ized" Bloch wavefunction for a quasi-particle should be 
taken as the eigenvector of spectral function. ^^'^■^ In this 
article, we provide a quite general derivation of the effec- 
tive EOM for quasi-particles based on the Keldysh for- 
malism, which confirms in part that this notion is valid 
in an arbitrary U{1) Fermi liquid. 

For N = 1 case, the form of this EOM is identical to 
Eq.([T]). However, in passing from the non- interacting ver- 
sion of gauge field definition, i.e. Eq.®, to that of the 
many-body case, we encounter an additional complica- 
tion. Because of the energy dependence of the self-energy, 
the renormalized Bloch wavefunction for quasi-particles 
- defined from the Green's function - depends on the en- 
ergy or frequency lo, in addition to crystal momentum 
k, i.e. \u°'{k,uj)). Accordingly, we are naturally led to 
introduce a sort of dual version of the electric field and 
associated electrostatic potential: 
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which, as well as the magnetic component, generally have 
non-trivial structures in the [d + l)-dimensional dual 
space. 

As will be shown in this paper, the renormalized k- 
space Lorentz force appearing in the effective EOM for 
quasi-particles is composed both of these magnetic and 
electric component estimated on shell: 
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Ect.fc is renormalized energy dispersion for a quasi-particle 
in question, which is defined as the pole of the single- 
point Green function with respect to u) (see Eq. (|23|) ). In 
this sense, this newly introduced electric component pro- 
vides another source of the Lorentz force in the k space, 
particular only to interacting Fermi systems. 

From the viewpoint of this EOM alone, however, the 
significance of the temporal component of the Berry's 
curvature is not clear. Namely, the renormalized Lorentz 
field can be also viewed as just the "magnetic compo- 
nent" defined in the codimensional space associated with 
the quasi-particle energy dispersion to — ea.k] 
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(see FigH]). Therefore, it is tempting to regard that the 
separate definition of the electric and magnetic fields, 
i.e. eqs. (|2l5p . is a sort of redundancy. Physically, this 
is somewhat sensible, since "sharply-peaked" quasiparti- 
cles ("infinite" life time with a definite dispersion rela- 
tion uj = ta^k) should be linked with those wavefunctions 
which are defined only in this codimensional subspace 
of the Lo ~ k Euclidean space. Such thinking suggests 
that only derivatives along this subspace are meaningful, 
which thereby involve appropriate linear combinations of 
LD and k derivatives. 

Nevertheless, it can be shown that the dual electric 
and magnetic fields defined in the uj-k Euclidean space 
do have distinct physical meaning. Specifically, we will 
show that these two fields enter into the linear response 
of the spectral weight to the applied electro and magnetic 
fields respectively (see eqs. (|84l88p ^. This conclusion is 
obtained by solving the quasi-particle spectral function 
perturbatively with respect to the gradient expansion. 
Therein, we observe, in both C/(l) and SU{2) case, that 
these Berry's curvatures characterize the 1st order correc- 
tion to the weight of the quasi-particle spectral function. 
Since the spectral weight in principle can be detected 
in a momentum resolved way, one may at least envision 
several photoemission experiments as candidate tools to 
make a "contour map" of the dual electro and magnetic 
fields separately (see sec. IV). 
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FIG. 1: uj-k Euclidean space and its codimensional subspace 
associated with u) = ea,k (Red line). B" can be viewed as 
just the "magnetic component" of Berry's curvature defined 
in latter subspace. 



A momentum-resolved measurement of the dual elec- 
tromagnetic fields is of interest not only as a mat- 
ter of principle, but also potentially in diverse ex- 
perimentally active areas such as anomalous and spin 
Hall effects in metals and semiconductors, and uncon- 
ventional superconductivity in some ferromagnetic met- 
als. Recent experimental activities in metal o^^'^^i^^ and 
semiconductor o — have focused on novel observations of 
anomalous and spin Hall effects, leading to controversy 
over the origin of these effects. On the one hand are a 
variety of proposals of intrinsic effects, related to Berry 
phases and the above-defined dual electromagnetic fields 
in the clean band theory limit^iiS. However, such Hall 
effects could also be of an extrinsic origin, i.e. a result 
of spin-dependent scattering of electrons from the spin- 
orbit potential of impurity atoms^i. Currently, the dis- 
crimination between these scenarios is rather indirect and 
based primarily upon the comparison to diverse theoret- 
ical model calculations and approximations whose physi- 
cal applicability is difficult to judgei^. If one could some- 
how experimentally "visualize" the distribution of these 
dual electromagnetic fields in a momentum space, inde- 
pendently from Hall/galvanomagnetic measurements, ex- 
perimentalists could readily judge for themselves whether 
intrinsic or extrinsic contributions are dominant in their 
material sample, by comparing with their transport mea- 
surements. 

A second potential use for a measurement of dual elec- 
tromagnetic fields stems for a recent proposal by Shi 
and Niu^ of a new formulation of the many-body prob- 
lem based on the non-commutative quantum mechanics. 
They begin with a quantization of the EOM given in 
Eq.Q, leading to a non-zero commutator between posi- 
tion operators, [Rj,Rm] = iejmil3i{k)^ , in addition to 
the usual commutators of momenta and position, i.e. 
[kj, km] = iejmi^i{R) and [Rj, km] = iSjm- As is the case 
with the momentum in the presence of a real magnetic 
field b, one can also introduce the canonical position op- 



erator R'j = Rj —iAj{k) such that their commutators be- 
come free from the dual magnetic field; [R'j,R'm] = 0. At 
a price for this, however, the electron-electron interaction 
acquires an additional phase factor, which can transform 
a repulsive electron-electron interaction into an attrac- 
tive one. Observing this, Shi and Niu attributed to the 
fc-space Berry phase an unique origin of the supercon- 
ducting phases realized only within ferromagnetic met- 
als, such as in UGe2, ZrZn2 and URhGe. A detailed 
experimental information constraining B{k) as a function 
of k, in combination with a measurement of band dis- 
persions would clearly provide these effective theoretical 
approaches an "ab-initio" model Hamiltonian, only to aid 
a quantitative comparison with material physics. 



The structure of this paper is as follows. In section II, 
we introduce a quite general multiple-band continuum 
model and the Keldysh equation for the lesser (greater) 
Green function/spectral function. In section III, based 
on this (dissipationlcss) Keldysh equation, we provide 
a systematic procedure for carrying out the projection 
into the low-energy band, so as to obtain the effective 
(reduced) Keldysh equation for those Green functions of 
quasi-particles forming a Fermi surface. 



Thanks to the Fermi liquid assumption, this effective 
Keldysh equation can be further solved in favor of the 
spectral function, perturbatively in the gradient expan- 
sion. Then, we observe in section IV that the Berry's cur- 
vatures, i.e. both dual electric fields and magnetic fields, 
indeed enter into the 1st order correction to the renormal- 
ization factor for quasi-particles. Using this solution for 
the spectral function, we further derive in section V the 
effective Boltzmann equation for the occupation number 
of quasi-particles in the phase space. Based on this equa- 
tion, we can finally read off the U{1) effective EOM in 
interacting Fermi liquids. Section VI is devoted to the 
summary and discussions of the present paper. 



A number of appendices describe other topics useful 
in understanding the main text in more detail. For 
completeness, we briefly review the Keldysh formalism 
and our notations in the appendix A. The appendix B 
is devoted to demonstrating the logical consistency be- 
tween our derived U{1) effective EOM and the so-called 
Ishikawa-Matsuyawa-Haldane formula'^^ within the lin- 
ear response regime. Appendix C describes how a U{1) 
magnetic monopole "current" in lu — k space (which ex- 
tends the notion of magnetic monopole density in fc-space 
into the to-k space) determines the distribution of dual 
electromagnetic fields in the U{1) case. As a specific ex- 
ample which gives a quantitative idea of the significance 
of the many-body correction i.e. £, we also present in ap- 
pendix D some specific model calculations based on an 
electron-phonon coupling Hamiltonian. 
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II. KELDYSH FORMALISM 



A. effective continuum model 

We begin with a general semi-microscopic model in 
which the electronic spectrum is described by a fc • p type 
expansion about some (arbitrary) point in the Brillouin 
zone. One may keep as many bands as are deemed close 
enough in energy to be relevant to the physics, and our 
arguments do not depend upon the order in the expansion 
in k. Familiar examples would be the multiple-band Lut- 
tinger models, Dresselhaus models, and Rashba model 
commonly studied in semiconductors, in which the ex- 
pansion point is chosen at the T point. The advantage 
of this formulation is that we can Fourier transform, in 
a usual way, to a continuous real space coordinate r, i.e. 
k — > — iVr. The non- interacting Hamiltonian is thereby 
expressed in terms of the slowly varying "envelope fields" 
-0a (f): 

no = Y. j driji{r)[H^{-iVr.r,t)]^^,4,^,{r), (11) 

where the band index a^'^ runs from 1 to N},. Specific in- 
formations about lattice, orbital and spin-orbit couplings 
are encoded into the matrix structure of [ffo]- Following 
the standard literatures^^, we will employ a short-ranged 
electron-electron interaction potential; 



^1 = E 



V ' '{r^ ro) 



X <('^i)^L(r2)V'4(r2)V'a;(?'i),(12) 



though the form of our results does not depend in detail 
upon this. 



B. dissipationless Keldysh equation 

Apart from the perturbation theory at equilibrium, 
the Keldysh formalism is constituted by the lesser and 
greater Green functions; 

g>(l|l') ^ -z(0(l)^t(i'))^ 

.g<(l|l') ^ ^{^\l'm)), 

where the time evolution of "0(1) = tpaiiri,ti) is deter- 
mined by the interacting Hamiltonian introduced above; 

Because these Green functions are not the time-ordered 
ones, they acquire the translational symmetry in space 
and time at equilibrium. Furthermore, they are always 
Hermite matrices with respect to space and time coordi- 
nates and band indices. 



Putting its derivation aside (see the Appendix A), let 
us begin with the Keldysh equation for these lesser and 
greater Green functions; 



(13) 



The (anti-)commutator here is defined by the convolution 
with respect to time t, space r and band index a; 



dl^(l|l) -(7(111') 



± 



J dlC(l|l) -5(111'), 



with J dl = J2ai I dridti. We denoted the bare Green 

function as Go which is composed of the quadratic part 
of the Hamiltonian; 



G„-i(l|l') ^ 



t,)5\n-r,). 



"d" above represents the spatial dimension of our system. 

Electron-electron interaction effects, on the other 
hand, arc encoded in the self-energy such as E^^(l|l') 

and S^^^''(l|l') , which are self-consistently given in 
terms of lesser/greater Green functions (see eqs. ()A5t - 
IA7P in the Appendix A). The former self-energy is the 
Hartree-Fock part, which is temporally instantaneous, 
i.e. E"^(l|l') - 5{ti - tv). The latter one is usu- 
ally dubbed as the lesser (greater) coUisional self-energy, 
which is at least 2nd order in electron-electron interac- 
tions and thus temporally non-instantaneous. The phys- 
ical role of this coUisional self-energy is two-folded, cor- 
responding to its following decomposition; 



a(l|l') 



r(i|i') 



1 tl-tl 



2 \ti^tv\ 



(S>(1|1')-S<(1|1')) 



i(Si^(l|l')-t-SA(l|l')), 

z(S>(l|l')-S<(l|l')) 
z(si^(l|l')-S^(l|l')), 



(14) 



(15) 



with a(l|l')* = a(l'|l) and r(l|l')* = r(l'|l). Ei^(l|l') 
and S^(l|l') in the right hand side denotes the retarded 
and advanced self-energy respectively. Thus, (t(1|1') 
stands for the real (Hermitian) part of the self-energy 
associated with the time-ordered Green function, which 
plays role of the renormalization of the quasi-particle en- 
ergy and wavefunction. On the other hand, r(l|l') is its 
imaginary (anti-hermitian) part, bringing about a finite 
life-time for quasi-particles. 

In a same way, we can define the real/imaginary part 
of the Green function; 



b(i|i') ^ 



■(l|l')-g<(l|l')) 
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^ ^(g«(l|l') + g^(l|l')), 
A(l|l') ^ *(g>(l|l')-g<(l|l')) 
^ ^(g«(l|l')-g^(l|l')), 



(16) 



(17) 



with b(l|l')* = b(l'|l) and A(l|l')* = A(l'll). Note 
that especially the latter one is nothing but the spectral 
function. As is clear from eq. (fT5|) . the Keldysh equa- 
tion for this spectral function is composed only by the 
real/imaginary part of the self-energy and Green func- 
tions introduced above; 



[Go' 



^HF 



0. (18) 



We will dub Go^(l|l') - Y:^^{l\r) - cr(l|l') as a 
"Lagrangian" L(l|l') in a sense that, when Fourier- 
transformed with respect to its relative coordinate, i.e. 
1 — 1', it reduces to lo minus a "renormalized" Hamil- 
tonian for quasi-particles (q.p.). Namely, its eigenval- 
ues specify the renormalized energy dispersions for q.p. 
as their zeros with respect to w, while their correspond- 
ing eigenvectors constitute an orthogonal set. Latter of 
which therefore can be regarded as (a periodic part of) 
the renormalized Bloch wavefunction. Accordingly, we 
are led to define the Berry curvatures and gauge connec- 
tions in the dual space in terms of the unitary matrix 
diagonahzing this Lagrangian (see eqs. ((39)) . p8)) and 
PT|) respectively). As is shown in this paper, our ef- 
fective Boltzmann/Keldysh equation legally derived via 
the projection process claims that the Berry's curvatures 
thus introduced indeed govern the effective EOM for q.p. 
in interacting Fermi systems. 

On the other hand, the commutator between the anti- 
hermitian part of the self-energy F and the hermitian part 
of the Green function b introduces a finite life time of 
q.p., i.e. the broadening of the spectral functions. When 
it comes to the q.p. closed to a Fermi surface, however, 
the Fourier-transformed F as a function of lu becomes 
as small as 0((w — jif'^T'^). Namely, the 2nd expres- 
sion of eq. dictates that, when analytically contin- 
ued from a Matsubara Green function, F at equilibrium 
is composed of a delta function, which imposes the en- 
ergy conservation on its internal lines. As a result, the 
momentum integral regions associated with the internal 
lines are restricted to be only near a Fermi surface at low 
temperature^. 

F in those electron-boson coupled systems with bosons 
having a finite excitation energy gets even smaller than 
this power-low decay. When jw — /x| is much smaller than 
this excitation energy wq, (the lowest order) perturba- 
tive calculations readily show that F vanishes exponen- 
tially, such as e"!'^"''!/'^" or e""^/"" (see Appendix D). 
Anyway, in both of these two cases, the broadening of 
the q.p. spectral function at sufficiently low T is by far 
smaller than the thermal line-broadening of the spectral 
function (~ T), which validates the so-called adiabatic 
assumption of Fermi liquid theory. 

Meanwhile, the hermitian part of the coUisional self- 
energy a remains finite even on a Fermi surface (w = /i) 



at zero temperature (T — 0). Namely, eq. indi- 
cates that, when analytically continued from a Matsub- 
ara Green function, a at equilibrium is composed of the 
principal integral rather than the delta function and thus 
free from the energy conservation imposed on its internal 
lines (see the appendix D for some example). As such, 
apart from the life time part F, its momentum integrals 
region for the internal lines are not restricted near the 
Fermi surface, which even causes the ultra-violet cut-off 
dependence of ct at T = 0. 

Because of these two different features generic in F 
and CT, we ignore in this paper the (intrinsic) lifetime 
effect F, while fully take into account the renormalization 
effects due to a. Namely, instead of eq. (fT5|) . we begin 
with the following dissipationless Keldysh equation for 
spectral functions: 



0. 



(19) 



Even this SU{Nb) dissipationless Keldysh equation is 
still non-trivial, due to the presence of the hand index, 
whose effect is the central issue of this paper. In the fol- 
lowing, we will present a general method of projecting 
out irrelevant band degrees of freedom associated with 
fully occupied bands and fully empty bands, so as to de- 
rive perturbatively the reduced Keldysh (kinetic) equa- 
tion only for the relevant bands constituting a Fermi sur- 
face. 



III. REDUCED KELDYSH EQUATION 

A. gradient expansion 

We will derive these reduced Keldysh equations per- 
turbatively with respect to the gradient expansion. The 
coupling constant of this expansion is a dimensionless 
quantity which measures how much a system is disequi- 
librated. To define this expansion accurately, notice first 
that the lesser and greater Green function at equilibrium 
acquire the translational invariance in space and time co- 
ordinates; g<(>)(l|l') = gafai,(r-i - ri>,ti - ty). Being 
given in terms of these Green functions, the lesser /greater 
self-energy and Lagrangian also become translationally 
invariant at equilibrium; L(l|l') = Lq^qj^, (ri — ry^ti — 
til). As such, when Fourier-transformed, the convolution 
encoded in the dissipationless Keldysh equation reduces 
a simple product only with respect to band index; 



La/3 
Aq/3 



drdt e-'«''+"* 



La/3 
Aa/3 



(r,t). 



Accordingly, we have only to diagonalize the Lagrangian 
so that an arbitrary diagonal A in this eigenbasis satisfies 
eg. pop at equilibrium. We will regard this trivial limit 
as the non-perturbed case and take into account disequi- 
librations as perturbations. 
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When a system is disequihbrated, the lesser/greater 
Green functions generally depend on the center of mass 
coordinate in space and time, i.e. R = ''^1,''^' and T = 
ti+ti/ . 

2 ' 

g<(»(i|i') = g<(>U^i;^,7^)- 

However, as long as a system is not so far from its equilib- 
rium case, their dependences on R and T are slowly vary- 
ing in comparison with the lattice spacing and inverse of 
band width respectively. As such, we could quantify the 
"distance" from equilibrium by a dimensionless ratio be- 
tween this slowly varying length (time) scale and a lattice 
spacing (inverse of the band width). To extract the lat- 
ter length/time scale, we have only to Fourier-transform 
the relative coordinate, i.e. r and t, so that the crystal 
momentum q and energy (frequency) cu are introduced^i; 

g<{uj,q;T,R) ^ j drdt e~'''-+''^'g<{r,t;R,T), 

g>{oj,q-T,R) = I J drdte-''^'-+'^'g>ir,t;R,T). 

Then, the derivatives of these Green functions with re- 



spect to q and Lu are quantities of the order of the Fermi 
length and inverse of a band width respectively. Thus, 
the (inner) products between the derivatives of these 
Green functions with respect to Q = {to, q) and those 
with respect to X = (T, R) can be regarded as a small 
dimensionless quantity, as far as a system is only weakly 
disequihbrated; 

r • g' » r ■ dgt » dxdx'g"" ■ dQdQ,g\ 

The superscripts "a" and specify the lesser or greater 
green functions. Since Lagrangian and spectral function 
are given in terms of these functions self-consistently, we 
can readily adopt the following relations also: 

L-A > dxl-dqk, dql-dxk 
> dxdx'\--dQdQ,A,---. 

Observing this, we expand the convolution in the (dis- 
sipationless) Keldysh equation, in powers of dqdx = 



[^x,^xX^Q,^Q,^]_ - {Xu ^ Qk]) - {X, ^ Q,] 
I 



(20) 



While we kept up to the 2nd order, one could explicitly 
write down the higher than this, by using the following 
formula for the Moyal product^^; 



{A®B){Q-X) 



d{l ~ 1') 



'•(1-1') 



J diAm B{i\i'), 



A{Q;X) B{Q;X) 



where d^ dg Og 



d^.{AB) 



x,v— ; = idx,dQ^A){dQ^dx,B). 
Note that the j and fc-summation in ea. (|20p run from 
to d, which will be made implicit from now on. The 
(anti-)commutators in eq. (|20p are taken only with respect 
to band indices, while the non-local correlation effect en- 
coded into the space-time convolution of eq. ()19|) is now 
perturbatively taken into account via the gradient expan- 
sion. 



B. projection process 

Our projection process is nothing but to solve eq. (PO)) . 
perturbatively in the gradient expansion, using the sub- 
stitution method. To be more specific, we are looking for 



A which satisfies this equation at a given order accuracy 
in the gradient expansion. To do this in a well-controlled 
fashion, we will first solve the off-diagonal elements of 
the spectral function in favor of its diagonal elements, by 
looking into the off-diagonal components of the matrix- 
formed KE given in eq. ((20|) . Substituting these solu- 
tions back into the diagonal components of the same KE, 
we therefore obtain sort of differential equations given 
only for the diagonal elements of the spectral function 
(sec. Ill B,C,D). Then, we will further determine appro- 
priate form of these diagonal elements, such that these 
differential equations are satisfied (sec. IV). 

As in standard perturbation theories, we begin with 
diagonalizing the zero-th order part. Introduce the uni- 
tary matrix which diagonalizes the Lagrangian L in the 
left hand side of eq. (^0)) : 



u^i U. 



(21) 



Then, a spectral function in this basis has only to be 
diagonal so as to satisfy eq. ipiH) at the zero-th order; 



A EE U^A U = 



Ai 







A 



(22) 
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We used "SansSerif for the spectral/green functfon rep- 
resented in the old basis, while "Roman" for those in the 
new basis. 

While each diagonal elements can be arbitrary at this 
level, from the physical point of view, they should be 
delta functions; 

A„^5(iJi) = Z(«)5(.;-e(°)). (23) 

Namely, e^^q above denotes the (rcnormalized) energy 
dispersion for the a-th band q.p., while Za stands for 
this q.p. spectral weight. The superscript (0) simply 
represents that they are quantities of the zcro-th order 
in gradient expansion. When employing this form as the 
0-th order solution, we will actually be able to satisfy 
the Keldysh equation up to the 1st order in the gradient 
expansion (see eq. (j74p ). Thus, we can justify posteriori 
that eq. is the appropriate 0-th order solution, based 
on which its higher order correction can be built up (see 
also the arguments in section IV). 

Under this unitary transformation, the usual deriva- 
tive encoded in Eq. (^0)) is replaced by the "covariant" 
derivative; 

U\dx^)U ^[Dx^B] = dxI3 + AxB-]3Ax-{'2'^) 

with B = ■ B • [/ and Ax = U^dxU- Namely, in terms 
of this derivative, our matrix-formed differential equation 
reads 

-[UA]='- [[bx, , id], [DqM ^ + • • • s F{A). (25) 

To find the spectral function satisfying Eq. ([25|) up 
to higher order in the gradient expansion, let us next 
look into the off-diagonal components of this covariant 
differential equation; 

— Ld^aAap + AapLdJj — Fap{{A^] , {A^^}). (26) 

Fap is a functional of a set of diagonal elements of the 
spectral function, i.e. {A^ = A^^^, and a set of its off- 
diagonal elements {^7,,} (7 ^ rj). Notice first that the 
right hand side of eq. ([^ is at least 1st order in gra- 
dient expansion. On the one hand, its left hand side is 
basically proportional to a direct band gap between a-th 
band and /3-th band, which is finite even at equilibrium, 
i.e. Ld^a 7^ Ld,(3- Thus the off-diagonal elements of the 
spectral function are of the order of 0{\dxdQ\), while its 
diagonal elements remain finite even at equilibrium. 

As such, we first solve these off-diagonal elements in 
favor of the diagonal elements of the spectral function, 
iteratively in gradient expansion. Specifically, to the 1st 
order's accuracy, we have only to replace the off-diagonal 
elements in the right hand side of eq. (I26p by zero; 

Aa^p = F^p{{A^}, {0}) • (Ldjj - Ld,^y^ = A^^l. (27) 

Using this, one could further obtain the solution of Aap 
to the 2nd order accuracy; 

A^'li{A,}) = F^pi{A,},{Ai^'^}), (28) 





a-Hl 


01 




F . S . 




a 

a-1 












Aa_i(m) 





FIG. 2: A schematic picture of the spectral function for the 
a-th band and its neighboring bands. 

or higher than that, 

=^«M{A},{4';)}). (29) 

When n being infinity, this clearly becomes an exact re- 
lation between diagonal elements and off-diagonal ele- 
ments. 

Substituting these solutions into the diagonal compo- 
nents of the covariant differential equation, we then have 
Nb decoupled equations which are given in terms only of 
a set of diagonal elements of the spectral function; 

= F^{{A^},{A^^}) 

= F„({A,},{A,,^4-)({A^})}), (30) 

with Fa = Faa- Then, a remaining task is to deter- 
mine a set of Nb diagonal elements of A, such that they 
observe these equations. {A^} thus obtained in combina- 
tion with A^ri = ^T^''({^7}) 3.re in principle equivalent 
to the exact solution of the an original dissipationless 
Keldysh equation, i.e. eq. (PU)) . 

However, the exact solution is apparently impossible, 
since it would require us, for example, to perform the 
iteration of eq. (|29p at an infinite time. Thereby, we are 
going to indulge ourselves in executing this sequence of 
the process, up to a given order in gradient expansion. 
The highest order up to which we have succeeded in ob- 
taining eq. (|30p is currently the 2nd order. Up to this 
order, we can readily ignore the 3rd order gradient ex- 
pansion term (and higher than that) denoted by "• ■ •" in 
eq. (PO)) . Furthermore, Faa already containing at least 
one pair of dxdq, the difference between eq. ([771) and 
eq. (I28p ends up with the 3rd order contributions, when 
substituted into Faa as in eq. Thus, we are ready 

to use eq. (P7|. 

The final simplification we will employ is that, for 
given fc, R and T, each quasi-particle bands have spec- 
tral weights at energetically well separated regions from 
one another. To be concrete, let us refer to a band which 
contains a Fermi surface as the a-th band. Then, the di- 
agonal elements of the spectral function corresponding to 
the other bands, i.e. A^^a, have negligible weights at the 
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low energy region, i.e. \uj — ea\ — \uj — ^ min^ I^q/jI- 
This can be seen precisely at equilibrium, where A-y^a 
is sharply peaked at an energy region separated from fi 
by the direct band gap A^q, (see Fig. [5]). Even off equi- 
librium, higher order gradient expansion corrections to 
Aj turn out not to make additional incoherent weights 
other than the delta function we originally have at equi- 
librium. Namely, the corrections to A^ appear in eq. (|23p 
only as the energy dispersion shift and the additive spec- 
tral weight (see for example eq. ([55)) ): 



,(0)+g(l) 



Accordingly, as far as this energy shift, i.e. e^^^ does not 
change the relative position of each quasi-particle energy 
dispersions, we are posteriori allowed to replace A^^^a by 
zero, for |a; — ea| <C min^Ao,^. 

To summarize the simplifications possible at the 2nd 
order analysis, we have only to derive the following equa- 
tions; 



= Fa,{{Aa,, Ay^a. E 
^7*7 ~ F-yjj[{Aa, Aj^a 



0},{AW}), 
eO},{0})-(L,,, 



(31) 



In the next next subsection, we will substitute the latter 
into the former, so as to obtain the (differential) equation 
only for Aa. We will dub the equation thus obtained as 
a reduced Keldysh equation (RKE). In the section IV, we 
will further find a Aa satisfying this reduced Keldysh 
equation. Thereby, we actually observe that Aa thus 
obtained is sharply peaked at w = e^, while having no 
incoherent weights at high energy sides. This observation 
will support posteriori the logical consistency built in our 
prescription described above. 



C. SU(2) FLs and additional coupling constant 

In the argument of the previous subsection, we have 
implicitly assumed that a Fermi surface in question is 
composed only by a single band (especially eqs. (|21l22p ). 
In general, this is true either in those metals without any 
centrosymmetric lattice point or in ferromagnetic metals. 
In such FLs, there remains only a charge degree of free- 
dom, while the (pseudo-)spin degree of freedom at each 
fc-point is usually quenched. This charge degree of free- 
dom is then described by the spectral function Aq,, which 
is a scoZar quantity. Accordingly, the RKE, i.e. eq. (pij) . 
becomes just a differential equation for this scalar func- 
tion, which we can name as a U{\) RKE. 

On the one hand, in usual paramagnetic metals having 
a centrosymmetric lattice point, each /c-point is (at least) 
doubly degenerate, describing spin-degree of freedom for 
quasi-particles. Namely, eigenvalues of our Lagrangian 
are always two-folded at equilibrium. Thus, a Fermi 
surface is composed by two degenerate bands, which we 



could name as SU{2) FL; 



Ld,i 







Here 2 by 2 matrices Ld^-^ should be proportional to a 
unit matrix at equilibrium. 

As such, the spectral function which satisfies eq. ([^0]) 
at the zero-th order in gradient expansion has only to be 
block-diagonalized. 



A = U^k U = 



Ai 







A 



(32) 



with arbitrary 2x2 matrices A^ (7 = 1, 

Corresponding to this generic degeneracy at the zero- 
th order, we will derive the RKEs in SU(2) FLs in favor of 
these 2 by 2 spectral functions, i.e. (7 = 1, • • • , Nh), 
so that spin and charge degrees of freedom are treated 
on an equal footing. As will be shown later, the RKE 
thus derived becomes a 2 by 2 matrix-formed differential 
equation, which we will dub as SU{2) RKE. 

Apart from relatively minor modifications, the deriva- 
tion of the SU{2) RKE also goes along with the same 
procedure as in U{1) case. Specifically, it also begins 
with the inter-band components of KE, which now take 
a 2 by 2 matrix-form; 



A 



J7]\ 



{{A,}, {A-,,,}), (33) 



where Aj^i is a 2 by 2 matrix, connecting the 7-th band 
and 77-th band. 

Being proportional to a unit matrix at equilibrium, 
2 by 2 matrices Ld.-y may be decoupled into its zero-th 
order part and a small degeneracy lifting part; 



(34) 



A finite ij = e-ycfz is generally originated from weak R- 
and T-dependences of the Green functions. Thus (di- 
vided by a characteristic band width) should be treated 
as a same order of quantity as Idxdg]. To convince our- 
selves of this more directly, consider a specific situation 
in the presence of a small magnetic field b, or equiv- 
alently a slowly varying magnetic gauge potential. In 
such a case, the Zeeman coupling energy between (bare) 
spin and b clearly should be included in e^. On the one 
hand, |9x9q| contribution turns out to be proportional 
to a spatial derivative of the external gauge field, which 
is therefore proportional to b also (for example compare 
eq. (|15)l with eqs. (|70ll7ip V As is obvious from this ex- 
ample, Idxdgl = Ai and e/A = A2 should be treated as 
same order of quantities; 



Ai ~ A2 



A. 
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Observing the inter-band components of the covariant 
differential equations, i.e. eq. we first relate the 

inter-band elements of A with a set of Ni, intra-band 
elements of A. To the 1st order in Ai or A2, i.e. C(Ai, A2), 
we have 



F^,{{A^}AO}) 



id) 



Substituting these 1st order solutions into the intra-band 
components of KE, we then obtain the following SU{2) 
RKE to the accuracy of ©(A^, A|, A1A2); 



\_Ld,a, Aaj 



I 



(0)iN-l 



0})-(i,,,-L^^!l) 



(35) 
(36) 



at equiliburium 



of f-equiliburium 




a-l 



FIG. 3: A schematic picture of the energy dispersion in 5*17(2) 
FLs. (Right): Each dispersion is doubly degenerate at equi- 
librium. (Left): When a system is weakly disequilibrated, 
the doubly degeneracy at each fc-point is lifted. The associ- 
ated splitting energy, i.e. e, is however much smaller than the 
typical band gap/width at equilibrium, i.e. A. 



The intra-band elements of A for the bands other than 
the a-th band, i.e. A^^a, were already replaced by zero, 
because of the same reason as we argued in the C/(l) case. 

In the next subsection, we will calculate eq. ([55]) in 
combination with eq. (|36p more explicitly, so as to obtain 
an actual form of the SU{2) RKE up to the order of 
©(A^). Out of the SU{2) RKE thus derived, one can 
immediately obtain the RKE in U{1) FLs, by regarding 
Aa as a scalar function and putting to be zero. 



D. Actual Derivations of SU{2) RKEs 

In this subsection, we will perform the actual calcu- 
lation of eqs. ([55)) and in a covariant way, only to 
arrive at the SU (2) RKE associated with the a-th band 
in question. By referring "in a covariant way" , we mean 
that every step in the following manipulation does not 
change its explicit form either under an SU (2) rotation 
Va within the a-th band or, under an SU{2Ni, — 2) rota- 
tion Va within its complementary space; 



U 



U ■ 



Va 

6 % 



(37) 



This is because any steps of the actual calculation are 
composed either by the inter-hand covariant derivatives 
or by the intra-band covariant derivatives^ both of which 
will be accurately defined in the next subsubsection. In 
addition to their definitions, some formula frequently 
used in the next next subsubsection will be also sum- 
marized in advance for clarity. 



1. arithmetic preliminaries 

We have already defined the SU{2Nb) gauge field and 
associated "covariant" derivative in eq. ([M]) . In an anal- 
ogous way, the SU{2) gauge fields and covariant deriva- 
tives for the two-fold degenerate a-th band are defined 
as follows; 



(38) 



where X could be any coordinates in the phase space. We 
sub(super)scribe a or a such as Ba or Bg,, only to suggest 
that this matrix is a 2 by 2 matrix associated with the 
a-th band or a {2Ni, — 2) x {2Nt, — 2) matrix in its com- 
plementary space respectively. When sub (super) scribed 

"aa" such as B^^a^ j ^ matrix ought to be regarded as a 
2 X {2Nb - 2) matrix. 

Apart from the derivative in the full Hilbert space, the 
covariant derivatives in its subspace do not commute with 
one another in general, 

= tdxA3c- - idx'A^, + ^ [i^ , A3c-] ■ (39) 

As is clear from the latter expression, (ixx' would be 
identical to zero, if this doubly degenerate a-th band were 



to subtend a complete set, i.e. X^o-^d 



1. 



Note that this derivative and the associated curvature 
^xx' transform in a covariant way under the unitary 



10 



transformation in eq. ([57]) : 

[D3c,Bo.] - vi[D'^,B^]vc., (40) 
while the gauge field is not; 

•^X ^i-^X^a + vidxVa- 

In addition to this gauge field in the a-th band, we can 
also define the inter- band gauge fields as the 2 x {2Ni, — 2) 
off-diagonal blocks of Ax', 



Ax = U^dxil = 



A\ A% 



X 



■^X ■^X 



Under the unitary transformation defined in eq. (|37p , the 
inter-band gauge field clearly transforms in a "covariant" 
way; 



^x 



while the intra-band gauge fields such as Ax and Ax are 
not. 

Using the latter ones, we can further define a 2 x {2Ni,— 
2) matrix-form derivative, which plays role of a sort of 
inter-band covariant derivative; 



dxBaa + Ax Baa ~ Baa Ax, 



where Baa stands for an arbitrary 2 x (2iV(, — 2) ma- 
trix. As long as this matrix is a covariant quantity, i.e. 
Baa — *■ vj^BaaVa, the intcr-baud derivative above clearly 
transforms in a covariant way under eq. (j37p : 



[Df,Ba^] ~*vi[DT,Ba6.]V&. 



(41) 



In terms of these quantities and derivatives, let us sum- 
marize henceforth several formula which become useful in 
the next subsubsection. Consider first a inter-band co- 
variant derivative of a inter-band gauge field; 



= [{dxU^) dx'U]^^^^^^,^ + [U'dj^x'tj]^^^^^^,^ 



(T" = ± 



+ E E [(dx'U^) U],a.^,.,,[u'9xU] 

J^a cr" — ± 

where (3 ^ a and cr, cr' = ±. Then, applying into the last 
two terms the following identity, 

E i"^>("^i = 1 - E E 

we can exchange the subscripts of the covariant derivative 
and the gauge field with each other; 

[DT,m,.\p.^)^[Dx%AT]i.\p.>y (42) 



In the actual calculations, this equality becomes very 
powerful, when combined with the "decomposition rule" 
of covariant derivatives such as; 

[Dx^BaaB'aB'l^J^ = [DJ" , Bas] -SgSg^ 

+Bac. [0% , Sa] B'L + Bac^K [DT,Ka] ■ (43) 

One should also note that curvatures either in the a-th 
band space or in its complementary space can be ex- 
pressed also in terms of inter-band gauge fields; 



-i{ATA%" -AS[fAT 
-i{ATAfi - A%^A 



T) 



(44) 
(45) 



2. SU{2) X SU{2Ni, — 2) covariant manipulations 

Using the arithmetics described so far, we will 
study eqs. (|35I36[) within the 2nd order accuracy in 
A. As is the case for a standard (such as Rayleigh- 
Schrodinger) perturbation theory, our 2nd order expres- 
sion for eqs. (|35l36p is given both by eigen-energies at the 
zero-th order and its eigen-wavefunctions (see eq. (j55p 
for example). Namely, differences between eigenvalues 
of L enter into a sort of "energy-denominator" , while a 
"numerator" in a usual perturbation theory is now tran- 
scribed into a gauge field (connection), which is noth- 
ing but the matrix element of our perturbation part (i.e. 
dxdg) among different eigenbases of L. Accordingly, just 
as in a usual perturbation theory, our 2nd order expres- 
sion for eqs. (|35l36p also depends on wavefunctions and 
eigenvalues not only for the a-th band in question, but 
also for the bands other than the a-th band. 

On the other hand, extensive semi-classical analyses 



in a non-interacting system 



5.6.7.8.9.10 



suggest that the 



low-energy effective theory for (quasi-)particles should be 
constituted only by Bloch wavefunctions and energy dis- 
persions for the a-th band in question, while free from 
details of the other bands. 

We shall show in this subsubsection that this is in- 
deed the case for the SU{2) RKE given in eqs. p5l36p at 
least up to 2nd order in A. To be more specific, we will 
transform eq. psp in combination with eq. (j36p into a 
more compact form rigorously up to 0{X^), only to find 
that they actually are given solely in terms of the SU (2) 

gauge covariant quantities such as i^xx' '^^'^ [Bxj ' ' ' ] • 
During this transformation, several formula described in 
the previous subsubsection, such as eqs. (|42ll45p . become 
very useful. 

To see this, let us begin with eq. ((35l) . i.e. the (a, a)-th 
component (diagonal component) of the original dissipa- 
tionless Keldysh equation; 

- [UAa]^Fi'\Aa) + Fi^H{Aa,}) + F('H{A,sm6) 

For the later clarity, the right hand side was decoupled 
with respect to different elements of A. This becomes 
possible clearly because our differential equation is at 
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most linear in A. Up to 1st order in Ai the first term L^; 
in the right hand side reads; — 



1 



Ld.a 

La,c 



(47) 



Observing eq. (|47p . notice first the commutator between 
where we have introduced following 2 by 2 hermite and Xi^ and the a-th band spectral function, i.e i„. This 
anti-hcrmite matrices respectively; 



commutator implies that the former hermitian matrix is 
^ the 1st order correction to the a-th band dispersion. No- 

Ma = -{Aq" {Ld,a — ~ i^j ^ Qj}],^^) tice also that N'a enters into the anti- commutator with 

Aa- As will be shown later, this anti-commutator lets 



(49) play a relevant role in determining the 1st order gra- 

dient expansion correction to the spectral weight (see the 
Ld,a denotes a {2Ni, — 2) x (2iVb — 2) diagonal block of section. IV for details). 



The 2nd order contributions in Fi^\Aa) and Fa'^'{Aaa) and F!f{Aa) are given as follows; 

■ { [Dq, J^q. , ^a] ] + ^gfiS: A^ + } 

[Dl, , Lfli - U.] - AT, [D%, , L^l - L,,^] - [D-x':, A^] {1^1 - 
.{ A-:] i. + A-: [D-Q^ , i„] + A-^ [D-Q^ , i.] } 
- {X, ^Q,\Xk,Qk} - {X,,Qj\Xk^Qk} + {Xj^Q,\Xk^Qk} - h.c. + ©(A^). (50) 



(2), 



- {X, ^Q,} - h.c. + 0(A2), 
2 



^'^S^^al^'lli - Ld,^)AT, - {Xj <-> Q,} - h.c. 0(A2 



(51) 
(52) 



where "• • •" in F'^^^ stands for those terms explicit in the right hand side of eq. (|T71) . To rewrite the latter two, i.e. 



' (2) " (3) 

Fa and Fa , only in terms of Aa, we first solve Aas and Aa in favor of Aa', 



AUAa) = ^{-AaA'^Q:{dx,LZi+[D%„L,,,]) + [D^^,Aa]ATdLd,.~LZi)}-{h,~LZi)-' 

- {Xk^Qk} + 0(A2), (53) 

A^s{Aa) = ^(Ld,,i:^"^iai5t--^r,^a-^Ql^rf,5)-(^<i,5-id,r,)-' " {Xk ^ Qk} + 0{X^), 

= -'-ATJcA^qI {Xk^Qk} + 0{\^), (54) 

with r\,b ^ a. One can obtain these relations, by looking into the inter-band components of the dissipationless 
Keldysh equation, precisely as in eq. ((36|) . We then substitute these two back into eqs. ()51|52p . only to obtain 
and Fa'' as a functional of Aa ; 

Fi'^ {Aa) ^ ^ LZ)AaA2^: [D%^ , L<Zi + L,,4 [L^l - L^.a] "^1^ + ^ {dx, L^) [^q, , Aa] A^A^ 



1 r 1 1 "ir^p^ r ^ (0) ^ I f ] 4Qa ^ A \ rS^^^K yjcKCS] f j (0) 

(0) 



4 , [D1, , + AT, - -Aa [Dt^^AfJ^ [D%^ , L^Z^ + L,,,] AT, 



aa. 
X, 



-(0) 
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—^A^^t [^1' + La..] [D%^. {lZ^ - L,..) ] {lZ^ ~ La,.)A% 

-{Xj ^ Qj\Xk,Qk) - {X,,Q,\Xk ^ Qk) + {X, ^ Qj\Xk ^ Qfe} -h.c. + ©(A^), 



(55) 



and 



taa Aaa a a. 



Qk y^d,a 



- [Xj ^ Qj\Xk,Qk] - {Xj,Qj\Xk ^ Qk] 

+ {X, ^ Qj\Xu ^ Qk] - h.c. + 0(A3). (56) 



Eqs. dnOl), (1S3 and in combination with eq. (gT]) 
are all the terms that enter into the right hand side of 
eq. dH]) up to 0(Af, AiA2,Ai). . 

Now that we have obtained an exphcit form of eq. , 
we shah next find out a 2 x 2 matrix Aa that satisfies this 
equation. Before doing this, however, it would be clearly 
helpful to simplify these more than 100 terms. In fact, 
we can further transform all these terms precisely into 
eqs. (|67ll69p . whose transparent expression helps us to 
find Aa in sec. IV. To be more specific, we can find, for 
any single term in eqs. (|50l55l56p . several counterpart 
terms with which it constitutes a certain SU (2) gauge 
covariant quantity enumerated in eq. (|67ll69p . 

Terms proportional to [ZJq,^^] To sec this ex- 
plicitly, focus first on those terms in eqs. (|50l55l56p that 
are linear in [_D", Aa\\ 

\[At, [D% , ~ LaA + At, [Dl, , L^l - U.] 
+ [Dt.AtMLZ'^-La..)] 

+ \{dx,LZ)[D2ik.MAZA'ir, 

- \ [Dq, , M A-qI [b%^ , + La A A% 

+••■. (57) 

where we have already used [i)Q,(L^*'^i — Ld.&)^^\ ■ 

iLf^~La,.) - -(4"ii-L,,s)-i-[I)a',(4"ii-i,.a)]. 
Note also that "• • •" above indicates those 3 kinds of coun- 
terpart terms with X and Q exchanged and their hermite 
conjugate terms; 

" • • • " EE -{Xj ^ g, \Xk, Qk} - {Xj,Q, \Xk ^ Qk} 
+{X^ ^ Qj \Xk ^ Qk} - h.c. (58) 

When taking these terms implicit into account, we may 
rewrite the 1st term in Ea. (|57p into the following; 

l{At [D% , L^i - L,,,] + AT, [D% , Lfli - L,,,] 



[L>x':,At]{LZ^l~L,A} 



- 4 [Dq, , M A-qI , L^i - La.,] A% 
\ATAD%.LZi-Ld.AA%l[D-QrM 



+ i [D%,ATA {LZ^ - U.)A%1 [D-Q^ , A„] + ■ • • . 

Namely, we replaced several terms in the left hand side 
by either their counterparts or their hermitian conjugates 
implicit in "• • •". We also used eq. (|42]) . only to obtain 
the 3rd term in the right hand side. Since we can regard 
that these terms were just swapped among their 3 other 
copies and their hermtian conjugates, we can begin with 
the following, instead of eq. ((57|) : 



Eq. ([571 
1 



^T;{dx,LZ) [D-Qk,Aa\ (iJtiS^ - igfit) 



-AtAD%.L,.ai-L,.a]A%l[D'^Q^,Aa] 



I [Dx%At] (4"ii - U.)A%1 , Aa] 



where "• ■ •" defined in eq. 

Compare the 2nd, 3rd and 4-th terms above with the 
decomposition rule for derivatives, i.e. eq. (|43p . Namely, 
we can unify these 3 into a single a-th band covariant 
derivative term; 



Eq. m = -^{dx,LZ)[D"Q,,Aa]^Q.x 



1 



I [D% , AZ {lZi - U.)A%-^ [D-Q, , Aa] 



■•(59) 



When combined explicitly with its counterpart with Xk 
and Qk exchanged, the 2nd term above is expressed by 
M.a (see eq. Accordingly, eq. ([FT)) can be rig- 

orously transformed into the following SU{2) covariant 
quantities, 

Eq. (EZl) = -'^{dx,LZ)[%,x,,[DQ,^M]^ 
~ {X, ^ Qj\Xk,Qk} - {X,,Q,\Xk ^ Qk} 
+ {Xj Qj\Xk ^ Qk} 

- ^[[D%,Ma],[D^,,Ao.]] -{Qj^X,}. (60) 
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Terms proportional to Aa Employing a similar 
manipulation, we can further simplify all the remain- 
ing terms, which are linear in either (a) Aa itself or (b) 
[Dq_, [Dq^, Act]]. As will be shown next, the latter one 
can be easily proved to be zero in total up to the order 
of 0{X'^). Thus, we will henceforth look into those terms 
in eqs. (|50l55l56p . which are proportional to Ac,; 



} 



— A A"" 



(0) 



- /I"" 4°" A /d"" ( T ^^"^ ^ — r J -] A"" 



(61) 



where "• • •" defined in eq. ([SS]) . We have already ne- 
glected those which are canceled either by their counter- 
parts or by their hermite conjugates. 

Notice first that, to 2nd order in A, we can regard that 
Aa in eq. I161\) commutes with other 2x2 matrices. This 
is because Aa reduces to a unit matrix at equilibrium, 
while eq. (pT|) is already at least 2nd order in Ai. Thus, 
commutators between Aa and other 2 by 2 matrices in- 
evitably end up with the 3rd order contributions in A, 
e.g. 

{dx,L<-^i)Aa[D"Q':,AT,]AQl 

= (ax.4°i) [DqI^ATMq"^- + 0(A?A2, A?).(62) 

Observing this, one can easily unify the 2nd term and 
(the hermitian conjugate of) 3rd term in eq. ([6T|) . so that 
they are given solely in terms of the Berry's curvature for 
the a-th band, 

1{AZ [D% , ijii - + {X, ^X,)} [D^q'^ , A^l] Aa 



2aa r j\a r (0) i 
(0) 



= ^{dx.LZ)ATM:^^Q"Mo.. 

Namely, when combined with its counterpart with Xj 
and Qj exchanged and their hermitian conjugate, 0{X^)- 
contribution of the right hand side above can be ex- 
pressed only by ^x,q/^ 



1 



" {^3 ^ Qj\Xk,Qk] 

= ^(5^.0[^5.-f^x,Q,]^a + 0(A^), (63) 



where we used eqs. (|62I43I44I) . 

The 1st, 4-th and 5-th terms in eq. (|6ip are described 
in terms of AAa and VLx q in total. To see this, let us 
begin with the 4th term; 



- A A°"^ 



= -l^aA-Q:{[D^^,[D%^,L^Z^ + U.]] 

- [D%, ^ [^S, ^ + Ld,4 ])aT,+ oix') 

= '^AaA-Ql [^kQ. > 4°ii - L,,4 + • • • 

= l^a [i^tit ' -^Q ■ (lTJ L,,,)A%^] _ 

~lAa [i^t-^S". AZ {L^Zi ~ U^)A%^] _ + 



•(64) 



where "■ • ■" stands for 3 other counterparts and their her- 
mitian conjugates. In the 1st equality, we have swapped 
1 /2 of the 1st term in the left hand side for its counterpart 
in (the hermitian conjugates of) {Xj ^ Qj\Xk ^ Qk}- 
In the 2nd equality, we have ignored 0(AiA2,Af)- con- 
tribution associated with the commutators between Aa 
and other 2x2 matrices. We further used eqs. (|39l45p . 
only to reach the final expression. 

Within C'(A^), the 2nd term in eq. is canceled by 
the 1st and 5th term in eq. 



(65) 



(1st) + (4th) + (5th) in eq. (HI]) 



1 1 



-Aa [A2,:AZ , A%^^ {L'Z^ ~ U.Mx^ 
-0(A=^) + • • • . 



Then notice that, when combined with its 3 other coun- 
terparts implicit in "• • •" , the coefficient of Aa above 
clearly reduces to the commutator between ^^g^Xfc ^"^^ 
M.a (see eqs. ()44l48p ). In combination with eq. (|63p. this 
dictates that C'(A^)-contributions in eq. (pT|) are indeed 
given only by AIq and 

eq. m = \{^x,^fM^Q.^^%Q^^Aa\^ 

- {Xk^Qk] ~\{{MaA%Q^_.Aa\^+0{\^). 

(66) 

Terms proportional to [Dq, , A^]] Those 
terms in eqs. (|50I55I56|) which are linear in the 2nd co- 
variant derivative vanish up to 0(}?\ One can eas- 
ily see this, by noting that any commutator between 
Aa and other matrices in eqs. (|50I55I56P ends up with 
0(A?A2,A?); 

- \ [dx,dx,Lfal^A%{Lf^\-U^,)AT, 
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1 



where "• • •" already defined in eq. 



3. 5(7(2) RKE up to the 2nd order accuracy and its 
physical implications 



According to the analyses up to the previous subsub- 
section such as eqs. (147160166^ . the C'(A^)-contributions 
of eqs. (|50l55l56p can be rigorously transformed into the 
following; 



(0)- 



(0)' 



(67) 
(68) 



(0) ■ 



(69) 



Ci above stands for the C'(A*)-contribution. Observing 
these results, notice that the first term in t2{Aa) 
is nothing but the anti-commutator between J\fa and 
Aa, encoded in 0(Ai)-contribution of F^^\Aa) (see 
eqs. (I47l49p ). On the one hand, the 2nd term in £2(^0) 
is obtained from the ©(Af )-contribution in F^^^ , F*^^^ and 
F'^', i.e. eq. (|66p . In spite of these apparently different 
origins, f^nd M^^, in these first two terms appear in 
a totally parallel fashion. We can see a similar feature 
also in the 5-th and 6-th term in £2(^0)- These observa- 
tions not only imply the consistency of our derived SU (2) 
RKE but also dictate that indeed plays role of the 
C'(Ai)-correction to the quasi-particle energy dispersion. 
Therefore, its non-trivial matrix structure as well as that 
of Cq, induces the spin-precession of quasi-particles. 

When a disequilibration is created only by the exter- 
nal electromagnetic fields, AAa defined in eq. (j48|) is com- 
posed by the spatial (magnetic) component and tempo- 
ral (electric) one. The former is the Zeeman coupling 
between an external magnetic field and internal mag- 
netic moment associated quasi-particle wavepacket. The 
later one is that between an applied electric field and 
internal electric dipole moment. To see these internal 
dipole moments explicitly, consider a situation in the 
presence of physical (i.e. external) electromagnetic gauge 
fields (ao,a), with corresponding physical electromag- 
netic fields b = V X a, e = Vijao — drSi- Then, fol- 
lowing the standard recipe, let us introduce the canoni- 



cal momentum k = q + a(T, R) and canonical frequency 
w' EE w-ao(i?)^; 

l"X;Q/ — l"fe=9+a,w'=w-ao/- 

Substituting this into eq. (|48p . one can then re-express 
Ma, such that it is given solely by the partial derivative 
with respect to these canonical quantities. Namely, this 
C'(Ai)-correction to a q.p. energy dispersion precisely 
reduces into an inner product between external electro- 
magnetic fields and sort of internal dipole moments; 

[Ma] = b • [M„] + e • [P„] 

x(%.4°ri4°i-L|5.,<)r'),(7o) 
-(a.„<ir|4°i-L|9.,4i"^')).(7i) 

Note that the magnetic dipole moment Mq, above 
reproduces eq. ([3]) in a non-interacting limit. On the 
other hand, the electric dipole moment has no non- 
interacting counterpart, since it is purely associated with 
the energy-derivative of the quasi-particle Bloch wave- 
function. 

As is trivial from its coupling with e, when considered 
in U{1) FLs, this electric dipole moment is literally 
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time-reversally even while parity odd; 

Pa(fc) = Pa(-fc) T-reversal, 
Pa{k) = — Pq(— fc) /— invserse. 

Therefore, in SU{2) FLs where both of these symmetries 
are guaranteed, becomes traceless just in a same way 
as Mq, does; 

Tr[P„(fc)] =Tr[M„(fc)] eeO. (72) 

When translated into the the SU{2) effective EOM, 
these two dipole moments, first of all, enter into that for 
the CP^ vector; 

= {Ma -b + P^ •e + ---}z, 

which describes the spin-precession due to the exter- 
nal electromagnetic fields. Observing eqs. (|67ll69p . note 
also that every derivative term found in a conventional 
Keldysh equation is now replaced by the corresponding 
SU{2) covariant derivative in our SU{2) RKE; 

Therein, its usual derivative part describes the charge 
degree of freedom, while the commutator with Ax stands 
for the precession of quasi-particle spin due to this gauge 
field. In this sense, our derived RKE treats the charge 
and spin degrees of freedom on a equal footing way, by 
using the SU{2) covariant derivative. 

To uncover the physical significance of remaining terms 
in £2 such as the 1st 4 terms and the final 4 terms, we 
further need to solve this reduced Keldysh equation in fa- 
vor of Aa (sec. IV) and to derive the effective Boltzmann 
equation (sec. V) respectively. 



However, it does not up to the 2nd order, i.e. £2(/oi) ^ 
0. To resolve this, we will introduce the 1st order correc- 
tion to the spectral functions; 

=/oi + /i + ---, (75) 

such that Aa satisfies the SU{2) RKEs even up to the 
2nd order accuracy; 

6 = £i(/i)-f £2(/oi)- (76) 

In general, we could further obtain the higher order cor- 
rection of the spectral function Aa, provided that the 
RKE is given up to the 3rd order's accuracy or higher 
than that, i.e. C = Ci -\- C2 + Cz + Ci A — ■. Namely, its 
3rd order part £3 determines the the 2nd order correction 
to Aa in terms of /i and /o; 

6 = A(/2)+£2(./i) + 4(/oi), (77) 

The following equations further specify the higher order 
corrections such as /a, • ■ •, in terms of /2,/i and /o iter- 
atively; 

6 = A(/3)+^2(/2)+4(/l)+A(/oi), (78) 

6 = 

However, our derived RKE being exact up to 
0{\\, A1A2, A2), the highest order to which accuracy we 
could determine A^ is /i. 

B. first order correction to Aa 



IV. PERTURBATIVE SOLUTION FOR RKES 

A. prescription 

When obtaining the reduced Keldysh equation de- 
scribed in the previous section, we have chosen the spec- 
tral function at equilibrium as a unit matrix with its 
coefficient to be a delta function of i',"-' ; 

Ac = <5(4"i)i - Zi'^Siu; 6("))i ^ /o(4"l)i- (73) 

This is not only because such solutions are physically 
sensible, but also because, starting from this zero-th or- 
der spectral functions, we can indeed keep on satisfying 
the RKE up to higher order in A. We will prove this 
point, by solving actually the derived RKE in favor for 
Aa perturbatively in A (see eq. ([ST]) '). 

Our choice of the spectral function at equilibrium 
clearly satisfies the derived SU{2) RKE up to the 1st 
order; 

= 6. (74) 



As we will show below, even up to this lowest order 
analysis, we obtain a non-trivial correction to the spec- 
tral function, having an interesting physical implication 
(next subsection) . Notice first the last 4 terms in eq. (p^ 
do not contribute at all to £2(/ol)j since they are anti- 
symmetrized with respect to the exchange between X 
and Q; 

A(/oi) = -y{[eo + ^o,^^x,Q,]_ + 
-^{dx^Lfl) i^^^.gj + ^{^Q,Lfl) [D'^x.A%q;\ } 

-/o{*(5Q,4°i) [Dx^etA - ^{dx,L^l) }. 

where Ct^a = Ca + ■M.a- Then, the last two terms above 
are readily canceled in eq. (fTG]) . when /i = —h,afo is 
substituted into Ci; 

A(-et,a/o)+^2(/oi) = 

+ *(5Q.4°i)[^l,f^x,Q,]}. (79) 
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Since et,a is regarded as a small quantity, —it.afo is noth- 
ing but the C'(A)-correction to the argument of /o, i.e. 
f(o) 

To set off also those terms in jC,2{fol) which are linear 
in /o itself, we need to introduce a correction to an over- 
all coefficient of /o, in addition to that to its argument. 
Namely, consider the following /i; 



1 



(80) 



Then, when entering into Ci, the 2nd term above totally 
cancel all the terms in the right hand side of eq. ((75|) . 

Notice that, this cancellation becomes possible, only 
because the relative ratio among the coefficients of first 
3 terms in C2{Aa) perfectly match with the ratio among 
the corresponding three terms in Ci{Aa)', 



1 - 

'\ . 
1 r 



4 

i 

+ 4 



Qk} + ---, 



^ £i(i„) = - [e„, i„] _ - [Mc,, i„] _ 

If either signs or coefEcients of any one of the 3 terms in 
C2 were not to meet with those of corresponding 3 terms 
in Ci, we could never have any /i satisfying eq. ([75)1 . 
In other words, there is no a priori guarantee for the 
existence of /i which satisfies eq. ([75)1 . In spite of this, 
both signs and coefEcients in C2 completely met with 
those in £1, we safely have eq. (l80l) as an appropriate 
C'(A)-correction to A^. Reversely speaking, this perfect 
coincidence between Ci and £2 indicates the validity and 
the consistency of our derived SU (2) RKE. 



C. Physical Implications 

Let us next argue the physical consequence of our so- 
lution Aa] 



Ar, 



/o(4°l)i - hMLZ) - ^^%Q,h{Lll). (81) 



(0) 



1 



(0) . 



As was already mentioned, the 2nd term above clearly 
stands for the renormalization of the energy dispersion. 
Namely, diagonalizing + A^a = we can transcribe 
this correction into that of the arguments of /o. 



/oi — ^t,af'o 



A finite €t,a,+ — ^t.a - thereby lifts the doubly degeneracy 
at equilibrium (Fig. HJ^b)). 



On the other hand, the last term in eq. ((8T|) changes 
the weight of the spectral function. To see this, let us 
integrate eq. ((8T|) with respect to frequency; 



(82) 



A^diO = X {1 - ^(f^X.Q,) 1^^,(0)}, 



where Za^ 



stands for the residue at L 



(0) 



0. This 



clearly dictates that, apart from the conventional renor- 
malization factor zj^'^ , the integrated spectral weight ac- 
quires an additional C'(A)-correction which is given by 
the SU (2) Berry's curvature. 

This curvature correction q ^ spin- selective ef- 
fect. Namely, choosing the basis diagonalizing this SU (2) 
Berry's curvature, we have. 



^ 2 ^ '■Xj Qj 





1 - h^x] 



(83) 



where the up-chirality electron and down one generally 
acquire the different correction to its spectral weight with 
each other. 

Observing eq. (jSSp . we can propose several photoemis- 
sion experiments as a candidate experimental tool to vi- 
sualize the dual SU{2) electromagnetic fields in a mo- 
mentum resolved way. To see this, consider the situa- 
tion in the presence of applied electromagnetic fields, i.e. 
b = V X a and e = Vao — OtSl. Then, using the canonical 
momentum and frequency introduced previously, eq. (I82p 
reduces into the vector product between the real electro- 
magnetic fields and the a sort of dual version of SU{2) 
electromagnetic fields; 

z„ ^ 4")x{i-l(f7j^.Q^.)|_^,o,} 

= X {l + ig".b + i£".e}. (84) 

Here the dual quantities in the right hand side are esti- 
mated on shell at equilibrium; 



£^^t{doA'^-d^A^+[A^,A'^]} 



(85) 
(86) 



(0),QCT|r, (0),aCT' 



[A^]^^,^{uZnd.u^^.l:-'), 



with m = 1, • • • d. 

As is clear from their coupling with real electromag- 
netic fields, B'^{k) and S^{k) as functions of k obey pre- 
cisely same symmetries as Mq, „i(fc) and Pa.m{k) do re- 
spectively. Thus, in a SU{2) FL having time-reversal and 
spatial inversion symmetry, these SU (2) matrices become 
both traceless; 



Tr[f«(fc)] ^Tr[g^^(fc)] ^0 



(87) 
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FIG. 4: A schematic picture of spectral functions in SU{2) 
case, (a): The correction due to the Berry's curvature changes 
its weight. Furthermore, the change for the up-spin weight 
and that for the down-spin are in general different, (b): Usual 
1st order correction change an energy spectrum. 



just as in eq. ([72]) . This suggests that the (inte- 
grated) spectral weights for the doubly degenerate quasi- 
particles, subjected under applied electromagnetic fields, 
acquire the 1st order corrections having a reverse sign 
with each other (see Fig. [3l|a)). 

The spectral weight for quasi-particle is in principle 
detectable in a momentum resolved way in terms of pho- 
toemission experiments such as angle resolved photoe- 
mission spectroscopy (ARPES) and resonating inelastic 
X-ray scattering (RIXS). The above theoretical observa- 
tion therefore leads us to raise the spin-resolved ARPES 
as the potential candidate tool to measure the dual SU (2) 
electromagnetic fields. Namely, observing the change of 
the spin-resolved spectral weight at each fc-point under 
small e or b, we can identify the dual electric or mag- 
netic field as its linear response. Furthermore, eq. ([57]) 
claims that, maximizing the measured linear responses 
with respect to the resolved- spin's direction, one can even 
determine the spin-quantization axis for the eigenbasis 
diagonafizing or f Thus, provided that these pro- 
cesses are performed for all the components of e and b 
separately, we could determine, as 2 by 2 matrices, all 
the 2 X d components of these SU (2) Berry's curvatures 
at an arbitrary k point on a Fermi surface. 

The (spin-resolved) ARPES experiment under exter- 
nal electromagnetic fields is, however, a difficult experi- 
ment. A more relatively less unrealistic proposal might 
be detecting abelian Berry's curvatures in U (1) FLs, such 
as ferromagnetic metals or paramagnetic metals without 
any centrosymmetric lattice point. In such U{1) FLs, 
we have only to regard all the SU{2) hermitian matrices 
in eqs. (|67ll69p as a scalar quantity, to obtain the U{\) 
RKEs and its solutions up to the 1st order accuracy^; 



1 



2 

= (l + i^"-b + f"-e)Z„^(cc;-e„). 

Note that the dual abelian electromagnetic fields are 
again defined to be on-shell as in SU{2) case; 

A-^^ = (89) 
with II = Lu,k. 



In such a J7(l) FL, the time-reversal (spatial inversion) 
counterpart being already lifted at equilibrium, these lin- 
ear responses against b and e can be measured only 
by the spectral weight itself for a given Fermi surface. 
Namely, the spin-filter are already implemented in mat- 
ters themselves. Thereby, in stead of ARPES experi- 
ments, the resonating inelastic X-ray scattering (RIXS) 
experiment, which is by far compatible with applied elec- 
tromagnetic fields, is a more promising photoemission ex- 
periment. 



V. U{1) EFFECTIVE BOLTZMANN EQUATION 
AND EOM FOR QUASI-PARTICLES 

In the previous section, we clarified the physical impli- 
cations of the first 4 terms of £2 given in eq. . We will 
now study on quasi-particle dynamics, only to interpret 
remaining 4 terms in C2 such as 



Co = 



(0)- 



{X, ^ Q,\Xk, Qk} - {X„ Q,\Xk ^ Qk} 
{X, ^Q,\Xk^Qk}. (90) 



Strictly speaking, when it comes to an EOM for quasi- 
particles, we in principle have to begin with the Keldysh 
equation for the lesser (or greater) Green functions g^^^-*, 
i.e. 



.<(»] 



(91) 



which is different from that for the spectral function A. 

As far as a system is not so far from its equilibra- 
tion, however, this difference could be regarded as small 
at T ~ and a; ~ /i. To be more specific, thanks to 
the boundary condition imposed on the Matsubara self- 
energy and its continuity, both lesser and greater coUi- 
sional self-energies '^^'-^{uj) above become small around 
w ~ (see the arguments in Appendix A2d and A3). As 
a result, we could begin with the following dissipationless 
equation, instead of eq. ([OT]) : 



(92) 



Note that a smaU S<'>(cj)|j^^^ does not necessarily lead 
to a small 'j(a;)|^^^. Namely, the latter one is given by 
the energy- (principal) integral of the former two; 



-\-oo 



dui' V 

271 UJ — L 



With a help of this biased treatment of the coUisional 
self-energies for different frequency regions, however, we 
can obtain clear physical interpretations for all the re- 
maining terms in C2 given in eq. (|90p . Namely, we can 
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apply precisely a same projection process onto eq. 
as we did for A. Thus, we begin with the exactly same 
reduced Keldysh equation for as in eqs. (|67ll69p . with 
being the (a, a)-th element of Nh times Nf, matrix 
g'^ = U^g"^ U . By constructions, this diagonal element 
g"^ can be decomposed into the product between the gen- 
eralized Fermi distribution Junction fa (scalar quantity) 
and the spectral function for the a-th band, later of which 
was already derived in the previous sectional; 



(93) 



Generally speaking, based on this decomposition, one 
must derive a coupled effective Boltzmann equation for 
spin and charge, out of the SU (2) RKEs for g^. However, 
since this analysis is somehow involved, we will hence- 
forth restrict ourselves to U{\) FLs, only to derive effec- 
tive Boltzmann equation for charge degree of freedom, 
i.e. EOM only for /„. 

In U[\) case, substituting eq. into the abehan 

RKE, we first obtain; 

= {Ci{Aa) + C2{Aa)] X /„ 

Notice that the spectral function A^ was determined such 
that Li{Aa) -)- C2{Aa) = (up to the 2nd order in gra- 
dient expansion). Thus, in the right hand side, we can 
safely drop those terms which are proportional to fa it- 
self. Notice also that the spectral function Aa thus de- 
termined is sharply peaked at w = . Thereby, we have 
only to integrate this equation over frequency, such that 
fa is replaced by the physical quasi-particle occupation 
number in the q-R space; 

na{q;R,T) = fa{q,uj = ea;R,T). 

The EOM thus obtained reads as follows; 

(1 - n^,^ + (dn.ea) n^,^ - (d.^ea) n%J drna 
= {dR^ea + (drea) nij,^ - (dR.ea) nin^ + n^n, 

+ {dq.Ca) ^R^rJ dq^Ha - {qj ^ 

which we can regard as the Boltzmann equation for the 
a-th band quasi-particles. Note that the partial q, R and 
T-derivatives encoded into the curvatures therein apply 
only onto their explicit dependences and do not apply to 
their arguments of e„, e.g. 



qiRj 



i^'qiRj)\i, 



Rj)\u 



As was explained in the introduction, however, one 
can also introduce the curvature defined in the co- 
dimensional space associated with uj = faiq^ R,T), e.g. 

Kr, ^ i{dqA%-dR^Al). 

= \u") = \u^)\^=,^. (94) 



Then the Boltzmann equation above can be also ex- 
pressed solely in terms of these curvatures in the q-R-T 
space. Specifically, normalizing the coefficient of drria, 
we obtain the following up to the 2nd order in A; 

= drUa + 

- {dq^e + n^^^ - {dR,e)ni^^ + idq,e)n%^JdR^na.{95) 

Accordingly, comparing this Boltzmann equation with 
the continuity equation, i.e. = 9t"-q + {dT<l)dq7ia + 
{dTR)dRna, one can readily read the effective EOM for 
the quasi-particle; 

^ = -dq^e - n^^^ + (dR^ea) ni^^ - (dq^ea) 



This effective EOM for quasi-particles proves, at least 
to the accuracy of 2nd order in A, that the EOM valid in 
"non-interacting" Fermi system^ can be also generalized 
into "interacting" Fermi systems with those curvatures 
defined in the co- dimensional space, such as eq. (j94l) . 
However, observing that the 1st order correction to the 
renormalization factor is characterized by the Berry's 
curvature in Euclidean q-u-R-T space instead of that in 
co-dimensional space (see eq. ([SS]) ). studying the 3rd or- 
der correction to this EOM is still an interesting open 
question (see also sec. IV). 

As an immediate application of this effective EOM, we 
can again consider the case with external electromagnetic 
fields; 



dR 

dT 
dk 

dT 



b X 



dk 

dT' 
dR 

dT' 



where — dk^a and reads. 



(96) 
(97) 

(98) 



Note that B" and are curvatures in Euclidean k- 
uj space (eas. (|2l3p and eas. (|5l6p respectively), while 
is the curvature defined in the co-dimensional space 
(eqs. (|9ll0p ). Observing eq. ([M)) . one then see that the 
intrinsic AHE in U{1) FLs should be defined in terms of 
both electromagnetic fields in the dual space; 



dk 

(27r)° 



B° 



X v„);nf(e„,fe),(99) 



where nf denotes a Fermi distribution function. 

Notice that the terminology "intrinsic" is now gener- 
alized into a slightly wider sense. Namely, the above ex- 
pression for the Hall conductivity contains not only the 
"fc-space magnetic field" effect, which is already present 
in a "non-interacting" limit, but also the "many-body 
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effect" via the corresponding "electric field" component. 
Based on the Ishikawa-Matsuyama formula and Fermi 
liquid assumptions, one can also see that this "many- 
body" correction indeed takes over (a part of) the so- 
called vertex correction to the static transverse conduc- 
tivity (see appendix B). This consistency check at the 
linear response regime strongly supports the validity of 
our derived effective EOM i.e. eqs.dM]) and (|97ll . 



VI. SUMMARY AND OPEN PROBLEMS 

To extract an information of low-energy effective the- 
ory in a generic multiple-band interacting Fermi systems, 
we have derived the Reduced Keldysh Equation (RKE) 
which effectively describes the charge and spin degrees 
of freedom for a specific single-band forming a Fermi 
surface. Our derivation is perturbative with respect to 
the gradient expansion whose coupling constant measures 
how a system is disequilibrated. It is, however, non- 
perturbative in the electron-electron interactions. In- 
stead, it relies only on the "adiabatic assumption" which 
is also utilized to validate the usual Fermi liquid frame- 
work. This assumption claims that, when it comes to 
the low-frequency region at sufficiently low temperature, 
the intrinsic life time of a quasi-particle (due to electron- 
electron interactions), which is 0{T^'^,{uj — /i)~^), be- 
comes much longer than the inverse of thermal broaden- 
ing of the spectral functions (~ 0{T~^) ). This assump- 
tion usually provides a finite energy region within which 
Fermi liquid theory works. Based on the same spirit and 
as a sort of zero-th order approximation, we completely 
ignored the life-time part (anti-hermitian part) of the col- 
lisional self-energies, while studying on its hermitian part 
on a general ground. 

Out of the RKE thus derived, we have succeeded in 
extracting several intriguing physical implications both 
in SU{2) and U{1) FLs. A first observation is that the 
linear response of the spectral weight for a quasi-particle 
with respect to an applied electromagnetic field, e/b, is 
characterized by (what we call) the dual electromagnetic 
fields, i.e. £°'/B'^. Since the linear response condition of 
small e and b precisely coincides with the condition of va- 
lidity of the gradient expansion, these linear response ex- 
pressions are asymptotically exact for any quasi-particle 
on a Fermi surface at zero temperature. Based on this 
theoretical observation, we also tried to give some rough 
idea of how to measure the SU (2) Berry's curvatures in 
terms of photoemission experiments, so as to provide 
some future directions for the (spin) galvanomagnetic 
community (see sec. IV). 

Another important achievement is our derivation of the 
effective equation of motion (EOM) for a quasi-particle 
in a U{1) FL up to the 2nd order accuracy in gradient 
expansion (see sec. V). From this equation, one sees that 
both the "fc-space magnetic field" and its correspond- 
ing temporal component, dubbed as "electric field", con- 
tribute to the Lorentz force which acts on quasi-particles 



in fc-space. In the linear response regime, this effective 
EOM is indeed consistent with the exact many-body for- 
mula for the static transverse conductivity (see appendix 
B). 

We conclude with a discussion of future problems and 
issues. First, one ambiguity remains unresolved in our 
derivation of the reduced EOM. Specifically, working only 
to 2nd order in the gradient expansion, we cannot exclude 
the possibility that in eq. may be more generally 
replaced by j^. These two forms are equivalent up to 
2nd order but not more generally. To fix this ambiguity, a 
further analysis including 3rd order effects in the gradient 
expansion would be necessary, which we leave for future 
work. 

A more pressing issue is the well-controlled treatment 
of the anti-hermitian part of the coUisional self-energies. 
Having completely ignored the corresponding terms by 
hand, our results are valid, in a strict sense, only for 
those quasi-particles exactly on a Fermi surface and at 
zero temperature. Since any physical quantities do in- 
volve non-zero excitation of quasiparticles (even in lin- 
ear response), this is clearly not satisfactory. While one 
might naively expect that coUisional effects can be in- 
cluded by direct analogy to classical kinetics (e.g. in the 
relaxation time or more sophisticated approximations), 
the possibility of some more interesting interplay between 
collisions and the Berry phase physics captured herein 
cannot be excluded. This may be of particular interest if 
one considers the instabihties of conventional metals to- 
ward ordered states such as superconductors, dielectrics, 
etc. A treatment of both the anti-hermitian part (decay 
of quasi-particles) and hermitian part (adiabatic trans- 
port of quasi-particles) of the RKE on an equal footing 
is certainly warranted in the future22i. 
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APPENDIX A: BRIEF REVIEW OF KELDYSH 
EQUATION AND COLLISIONAL 
SELF-ENERGIES 

To make this paper self-contained, we will briefly re- 
view the derivation of the Keldysh equationS^, defining 
several kinds of self-energies introduced in section II, such 
as S^^, a and S^'^- In addition to this, we will also 
look into the boundary condition especially imposed on 
the collisional self-energy E^''' . This boundary condition 
with additional arguments validates the dissipationless 
Keldysh equation for green functions, i.e. eqs. (|92IA26p . 
starting from which we derived the C/(l) effective Boltz- 
mann equation in sec.V. 

In the next subsection, we will begin with the time- 
ordered Green function defined on the imaginary time 
domain G(l,l') (see eq. (|A1[) for its definition). Specifi- 
cally, the Dyson equation for this temperature (Matsub- 
ara) Green function is derived first. The Hartree-Fock 
and collisional self-energy encoded there are also itera- 
tively defined in terms of this Matsubara Green function 
(see eqs. (|A5IIA7|) ). 

The time-ordered Green function as a function of imag- 
inary time t is analytic separately in two region: Im t G 
[—(3,0] and [0, /?]. In other words, we can introduce two 
functions, usually dubbed as lesser and greater Green 
functions, which coincide with this time-ordered one and 
which are analytic in these two regions respectively 



1 Tr[e-'3(«-f^)r{^ V(l)^^(l')}] 



G(l,l') = 



G>(1,1') 
-G<(1,1') 



for 
for 



Im ti < Im ty , 
Im ti' < Im ti, 



(compare eq. ([XT]) with eqs. (|A8IA9|) ). In the 2nd sub- 
section of this appendix, keeping the analyticity of these 
two, we will extend its time domain from the imaginary 
time into the real time domain. Correspondingly, the 
Dyson equation is also analytically continued onto the 
real time domains (see eqs. (|A20IA22IA23|) '). 

We will see in the final subsection that periodic bound- 
ary condition imposed along the imaginary time axis re- 
lates the lesser and greater collisional self-energies defined 
on the real time domain with each other, (see eq. (|A28p ). 
This relation guarantees that, at zero temperature and 
at equilibrium, lesser and greater self-energies have no 
weight at w < /X and co > fj, respectively, i.e. eqs. 
(IA24IA25|) . These observations allow us to approximate 
the Keldysh equation derived further into the so-called 
dissipationless Keldysh equation (see the arguments from 
eq. (|X22|) to eq. (|M6)) ). 



1. Dyson equation for temperature Green function 

Let us begin with the imaginary time Green function 
defined in the "interaction" representations: 



G(l, !';</), to) 



^ i (r{^^(i)^t(r)}) 

-^ {T{S}) 



S = exp 



Tr[e-/3(w-pA')r{^}] 
I I d2 0(2) -7^(2) . (Al) 



Note that the argument of the fermion operator "j" is 
an abbreviation of (rj, tj, aj). Accordingly, J^'j* dj in- 
cludes not only the integral with respect to the imaginary 
time tj , but also the summation over the band indices aj 
and the integral over the spatial coordinate; 



to-i/3 



to-i/3 



dti / dr 



"T" is the imaginary-time-ordering operator along 
[to, to — i/3] , where to is always real- valued (In the next 
subsection, we will set to to — oo, after analytically con- 
tinue ti and ti' onto the real time domain.); 



ipil) V(l') for Im ti < Im ty , 
-V'(l') V'(l) for Im ti, <lmti. 



The time dependence of the fermion operator ipil) 
ipaiiri,ti) is specified hy H = Ho + Hi: 

-idtM^) = [w,V(i)], 



with Ho and Hi given in eqs. (|llll2p . Then taking the 
ti-derivative of eq. (jAl|) . we have the following EOM for 
this 1-point Green function; 



'"di Goi(l,i)G(l,l') = 5(1-1') 

-I J d21/(l,2)G2(l,2,l',2+)|t,=t,_,|,|, (A2) 

where the inverse of a bare Green function G^^ (1, 1') and 
two-point Green function are defined as follows; 

Go '(1, 1') = [{tdt, - m)Sa,a,, - Ho]Sil - 1'), 
G.(l,2,l,2)= . 

"2-I-" in eq. (|A2p means that its temporal argument is 
chosen to be infinitesimally later than t2 along the imag- 
inary time axis; 2+ = (a2,?'2,t2 — i|e|)- 

The auxiliary scalar potential 0(2) entering into 
eq. (jAip as an "interaction" is utilized so as to describe 
the above 2-point Green function in terms of the self- 
energy S and 1-point Green's function G; 



to-i/9 



dl [Go\l, i) - S(l, 1)] G(T, 1') = 5{1 - l')(A3) 



Notice first that the functional derivative of the 1-point 
Green function with respect to 0(2) brings about the 2- 
point Green function; 

G2(l, 2, 1', 2+) = [G(2, 2+) - ^] G(l, l'). 
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Thus, substituting this into cq. (jA2p . we obtained a 
closed equation for the 1-point Green function; 



to-i/3 



dl G-i(l,I)G(I,l') = <5(l-l') 



i / d2 V{1,2) G(2,2+) 



|t2=ti-i|e| 



G(1,1').(A4) 



This self-consistent equation can be readily transcribed 
into that for the self-energy E. Matrix- multiplying 
eq. (|A4p by G^^ = Gq ^ — S, we first obtain the following; 



E(l, 1') = -*<5(1 - 1') / d2 y(l, 2)G(2, 2+)|4,=*,_,|,| 



|t2 = tl-i|£| 



Using G-i = Gg-i - E and °"g^,^,V ' ^ --^(1 - - 



■?G-^(l,l') _ 
50(2) 

2) in the right hand side above, the following iterative 
equation for the self-energy S is derived: 



E(l,l') = -i5{l-l') yd2y(l,2)G(2,2-H)|t,=t,_,|,| 

-f i^(ii-iiOF(l,l')G(l,l')|v=*,-,|.| 

/E(i,l') 



^ i I dl I d2 V^(1,2)G(1,1) 

J to 



(A5) 



Being proportional to d{ti — ty), the first 2 terms 
are temporally instantaneous, which we usually call the 
Hartree-Fock term; 

= -iS{l-V) yd2y(l,2)G(2,2+)|4,=,,_,|,| 
+ iSih - tv)V{l, l')G(l, l')|t,,=ti_^|.|.(A6) 



A successive iteration through the 3rd term of eq. (|A5|) 
leads to the self-energy with higher order in V, which is 
clearly temporally non-instantaneous. We dub this non- 
instantaneous self-energy as the collisional self-energy 
henceforth; 



E,(1,1') = E(1,1')-E"^(1,1'). 



(A7) 



Eq. (|A3p in combination with eq. (|A5P is the Dyson equa- 
tion for the imaginary-time 1-point Green function. 



2. Analytic continuations 

a. analytic continuation of and 

Real-time Green functions are analytically continued 
from the following lesser and greater Green functions de- 
fined on the imaginary time domain; 



G>il,l';cj),tn) = -i 
G<(l,l';0,io) = I 



I 



(A8) 
(A9) 



U{t,t') = T^exp [ — 1 (il0(l) •n(l)] } characterizes the 
temporal evolution along the imaginary time axis [^o , to — 
if3] and thus is not unitary operator. 

Note that, at equilibrium and with = 0, these two 
functions are analytic in Im {ti —ty) e [—(3, 0] and [0, f3] 
respectively, regardless of Re ti and Re ty. Then, as- 
suming that these analyticities also hold true in a weakly 
disequilibrated system with finite 0, we will introduce 
separately the lesser and greater Green functions defined 
on the real time domain as follows; 

' &aiay{n,ri' : Re ii,Re ti,]. 



limto^_oo G<(1, 1'; 0, to)\ Im ti=Im tj,-|-|e| = |e 



ga^a,,irl,rl> : Re ii,Re tv; 
= limt,^_ooG>(l,l';0,to) 



IfAlO) 



I Im ti + |c|=Iin ti/=Oi 

by taking e infinitesimally small. Choose the test scalar 
field 0(1) such that it vanishes at Re ti — » —oo. Then, 
^(^Oi ^0 ~ appearing in both the denominator and the 
numerator of eqs. (jA8|A9|l reduces to unit, when to 



—oo. Accordingly, the two functions in the right hand 
side are solely defined on the real domain; 



•(1,1';0) = 
:(1,1';0) = 



(All) 
(A12) 



with ^^(1) = U\-oo,ti)ip{l)U{-oo,ti). Now that 
U{—oo, t) is a unitary operator, we can construct a real- 
time-ordered Green function in terms of eqs. ()A11IIA12|) : 



;(1,1';0) = I (r{^4l) ^^(1')}) 

g>(l,l';0) for ti < tv, 
-g<(l,l';0) for tv <ti. 



(A13) 



b. analytic continuation of and 

The Dyson equation being composed also of the self- 
energy, let us next look into the analytic continuation 
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of the coUisional self-energy. The analytic property of 
the collisional self-energy can be obtained from that of 
Green functions at equilibrium. At equilibrium, both the 
1-point Green function G(l,l') and the 2-point Green 
function 62(1, 2, 1', 2 + i|e|)|f2=ti-i|£| are analytic in two 
regions, i.e. Im {ti — ty) G [— /3,0] and [0,/3], separately. 
Compare these analyticities with the relation among the 
self-energy, 1-point Green function and 2-point Green 
function; 



to -10 



d2 V{1,2) G2(l,2,l',2-f-i|e|)|t,=t^ 



to 



\t2=ti-i\t\ 



to-ip 



d2 E(l,2)-G(2,l'). 



(A14) 



Then one can readily see that I]c(l,l') is also analytic 
both at these two regions, but separately. Thus, we can 
introduce the lesser and greater self-energy, such that 
they are analytic in the following two regions respec- 
tively; 



Sc(l,l')^ 



S>(1,1') for Im h < Im ir, 
-i;<(l, 1') for Im tr < Im h. 



Exactly in a same way as we did for the 1-point Green 
function, we can then formally extend these two self- 
energies onto the real time domain separately ; 



' ai ,ai, (''I ' ^1' • 1 1 , Re ; (/)) 

= Mmtg^^oc S<(1, 1'; (p, io)| Im ti=Im ti,+\e\ = \e[ 

^^ai,ai,(''i'^i' '■ ii,Re ti';(/)) 

= limfQ^_oo S>(1, 1'; (j), to)\ Im ti + \e\=Im tj^,=0- 



(A15) 



c. analytic continuation of Dyson equation 

Now that both the self-energy and Green functions are 
analytically continued onto the real-time domain, we will 
derive the Dyson equation for these real-time functions. 
We begin with that for the imaginary time domain; 

(Go1-EHF).G<(1,1') 

dl I]c(l,l) G(l,l'). (A16) 

'to 

Decompose first the right hand side into three terms, such 
that each term is expressed solely in terms of lesser / 
greater Green function / collisional self-energy. Namely, 
depending on whether ti locates within [to, ti], [ii, ^i'] or 
[tvyto — the right hand side can be divided into the 
following three parts; 

dl S>(1,I)G<(I,1') + / dl S<(1,1)G<(1,1') 

I to Jti 

r ' dl I]<(1,I)G>(I,1'), (A17) 
It,, 



where the integral path is still along [to,^o — iP] (see 
Fig. [5]^a)). Then, we will deform this integral path so 
that both ti and tii reach the real-time axis. During this 
deformation, however, to keep the analyticity for all the 
Green functions and self-energies in the both right and 
left hand sides, we have to observe the following condi- 
tion; 

- (3 <lmty < Im ti < 0. 

Due to this condition, for given "target" Re ti and Re tii , 
our integral path is uniquely deformed into a "L" shaped 
path depicted in Fig. Ob-c); 



(G„-1-EHF).G<(1,1') 



di I]>(1,1)G<(1,1') 



to 
to-iO 



*i' 



dl I]<(1,T)G<(1,1') + 



ti 



dl I]<(1,T)G>(1,1') 



to~if3 



to-iO 



dl I]<(1,T)G>(1,1'). 

(A18) 



Note that the 1st three integral paths in the right hand 
side are all along the real time domain, while the last one 
is strictly along [to, to — i(3]. So as to eliminate the final 
term, we will take to to be —00. In this limit, G^ (1, 1'; 0) 
in the last term should vanish, since |Re i j — Re ti' | ^00. 
This procedure simultaneously completes the analytic 
continuations of all the Green functions / collisional self- 
energies encoded in this equation (see eqs. (jAlOIAlSI) ); 



(Go 1 - 



<(!,!') 

dll]>(l,l)g<(l,l') + r'dlS<(l,l)g<(l,l') 
Jti 

+ I dlS<(l,l)g>(l,l') 



(A19) 



In a similar way, we can easily obtain an equation of 
motion for the greater Green function; 



r dl [E>(i,T)-E<(i,i)] g<(i,i) 

^ — 00 

dl I]<(1,1) [g>(l,l')-g<(l,l')]. 



1-1 vHF\ > 



(1,1') 



(Go^-S^ 

dl[l]r(l,T)-S<(l,l)] g>(l,l) 



dlS>(l,i) [g>(l,l')-g<(l,l')]-(A20) 



At equilibrium, both Green function and self-energy 
become translationally invariant in space and time, e.g. 

g<(»(l,l';0^O) = gSl(^i-'-i',ti-ti'), 
E<(»(1,1';0^O) = ^<(»^^^^Jr^-rr,t,-tr). 
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to- 



Ret^ Re ti 



to-iP 



(a) 



(b) 



Reti' Reti 



to-ip 



(C) 



FIG. 5: (a): The integral path for eq. (|XT6)) . (b): Integral 
path for eq. (|A18|) in the case of Re ti < Re ti' . (c) That for 
Re ti > Re ti' . 



Thus, when Fourier-transformed with respect to these 
relative coordinates, the convolutions appearing in 
eq. (|A20[) would have reduced to a simple product, if 
the associated integral region were to range [— c», oo]. In 
reality, however, their temporal integral regions do not 
range from — oo to +00. As a result, when Fourier trans- 
formed, they end up with several principal integrals with 
respect to the frequency. The hermitian part of the col- 
lisional self-energy and Green function solely take over 
these principal integral parts. To be specific, introduce 
the following two functions; 



b(l,l') 



'i^[g>(l,l')-g<(l,l')], 

[I]>(1,1')-S<(1,1')]-(A21) 



2\h 
1 ti 



tv 
tv 



2 \ti-h 



Then we can formally rewrite the above equation in such 
a way that their convolutions with respect to time are 
always taken over [—00, 00]; 



(Go 1 - E"^) ^< 



(7 (g) g< + S< «) b 



--(Er®g<-S<®g>), (A22) 



1 



(S>®g<-E<®g>). (A23) 



where we used the following abbreviated notations; 
{A®B)(l,\')= / dl 1)5(1,1'). 



d. dissipationless Keldysh equation 



Eqs. (|A22IA23|) are what we call the Keldysh equa- 
tion, from which we readily obtain eq. p^ '^^. Notice 
that the re-expression in terms of a and b is not just 
for its formality, but is in fact an important step to 
approximate this Keldysh equation, based on physical 
arguments. As we see in the next subsection, at zero 
temperature and at equilibrium, the periodic boundary 
condition along the imaginary time domain ensures that, 
when Fourier-transformed with respect to their relative 



coordinates, the lesser/greater coUisional self-energy van- 
ishes for those lo greater/less than chemical potential /i 
respectively; 



E^(g,a;) =0 for w > /i, 



for LO < fi. 



(A24) 
(A25) 



Thereby, provided that a system is weakly disturbed from 
this equilibrated situation, we could still assume that 
these relations hold true for the coUisional self-energy, 
with /i generalized into a function depending on R and 
T; fi fj,{R,T). Observing further the continuity of 
these coUisional self-energies around w ~ /x, we could re- 
gard both of them to be sufficiently small quantities at 
these low frequency regions, i.e. ~ /i. 

The hermitian part of the self-energy ^{lu), however, 
is not necessarily a small quantity, even for lu c:^ fi. 
Namely, when Wigner-transformed, those lesser/greater 
self-energies for high frequency regions do contribute to 
(t(w) at the low frequency regional; 



a{q,R;uj,T) = 



2n u! — uj' 



j{f:< + t>){q,R;io\T). 



Due to this different behavior between (j{u!) and 

T,^^^\lu) around w ~ /i, we can approximate 
eqs. (IA22IA23[) into the following dissipationless Keldysh 
equations; 



(Go 1 - - a 



i<{>) 



0. 



(A26) 



which is supposed to be still valid for a; ~ /x at sufficiently 
low temperature. 



3. periodic boundary conditions 



Observing eq. (jA14p . we will first read the boundary 
condition (b.c.) for the coUisional self-energy, out of 
those for the 1-point and 2-point Green functions. No- 
tice that the 1-point Matsubara Green functions defined 
in eqs. (jA8IIA9l) obey the following boundary condition 
at equilibrium; 



G<(1,1') = ^ef'''G> {1-1/3,1') 



(A27) 



for — /3 < Im ty < Im < and arbitrary Re ti and Re 
ty . In a similar way, we can also see from its definition 
that the 2-point Matsubara Green function observes, 

G2(l,2,l',2-|-j|e|)|4,^4,_,|,| 
= -e^^G2{l - if3, 2, 1', 2 -I- i\e\)\t,=t,-^p-^\e\ 

for —(3 < Im ti' < Im ti < and arbitrary Re ti and 
Re ti'. When compared with eq. (|A14p . these two then 
lead to the following boundary condition for the coUi- 
sional self-energy; 

E,(l,l') = -e^''I]c(l-*/3,l')- 
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While being imposed on the imaginary time direction, 
this equation becomes useful, when Fourier-transformed 
with respect to Reii — Reti»^; 

I]<(ri,ri,;c.) =e-''(-^)E>(ri,ri,;^). (A28) 

Namely, eq. (jA28|) requires both of them to strictly ob- 
serve eqs. (IA24IA25[) at zero temperature, since either 
T,f or Sj* can not be singular in a usual metal. 



which is known as the Ishikawa-Matsuyama formulaSi*^. 
In this appendix, we will show that our derived expres- 
sion for the C/(l) Hall conductivity, i.e. eq. ((99)) . is in in- 
deed consistent with this many-body formula, using only 
Fermi liquid assumptions. Namely, we will assume that, 
when diagonalized, each eigenvalue for the time-ordered 
1-point Green function has a pole sufficiently closed to 
the real-axis, such that the corresponding quasi-particle 
life time is infinitely long; 



APPENDIX B: CONSISTENCY WITH 
ISHIKAWA-MATSUYAMA-HALDANE 
FORMULA 

The conductivities being given by the current-current 
correlation functions, they are usually subject under the 
vertex corrections in interacting Fermi systems. How- 
ever, when it comes to the static and transverse com- 
ponent, the Ward identity relates this vertex part with 
the derivative of the 1-point Green function, such that 
the Hall conductivity is expressed solely in terms of the 
single-particle Green function; 



2h {2TrY 



dk 



2tt 



-duj dkn dk,. 



(jj — Ea,ki^) — isign(a; — /i) • 0+ 



We can identify U above as the unitary matrix diagonal- 
izing our Lagrangian L. 

Specifically, let us choose the basis in eq. (|B1[) . such 
that the Green function is diagonalized; 



e exfj.1, 
2h (27r)'i 



dk I — e^-o^ 
2tt 



Tr[{(a^C/) ga U'^ + U (d^gd) U'^ + U ga {d^U-^)} 

U 9d U-'{idk.U) ga' + U {dkj^') + U g^' idk..U-')} 
I 



(B2) 



These 27 terms can be classified into 3 types, according 
to the matrix structure within the trace. To see this, no- 
tice first that, irrespective of being differentiated or not, 
(jd and g^^ alternate each other within the trace. When 
we have either [/-1(9qC/) = Aq or {dqU-'^) ■ {dq^U) be- 
tween a pair of neighboring gd and g^^, this pair of green 
fimction and its inverse can not be directly- connected. On 
the one hand, when both U and between a pair of 
gd and g'^^ are free from the derivative, these two can be 
clearly directly-connected. We will first classify all the 
terms appearing within the above integrand in terms of 
the number of directly-connected pairs of gd and g^^- 

One class is those terms having no pair of directly con- 
nected Green function and its inverse. For example, the 
following term belongs to this class; 

Tv[U{d^gd)i}-' Ugd ' {dk, U-')U9dU-\du^ij)9^'ij-'] 
= Tr[(a^g^i) Av 9d AJ, 

where one pair of gd and g'^^ was directly-connected only 
to reduce into a unit, i.e. gd ■ g^^ = 1. Such an anni- 



hilated pair is not regarded as a directly-connected pair. 
Among 27 terms enumerated above, there exist 16 terms 
belonging to this class, all of which can be summarized 
into several total derivatives; 

dk^[TT[A^ gd A, gd']} 
+ 9fc^{Tr[Av ffd Ag^']} 
+ a^{Tr[A„ gd A, ffd'll^O. 

Being contour-integrated, all these terms vanish as indi- 
cated. 

The 2nd class is those terms having two directly con- 
nected pairs. Among 27 terms enumerated above, there 
clearly exists only one such a term; 

Tr[id^gd) [dkJd') 9d {duJd')] 

= Tv[{d^gd) {duJd')9d {dkjd')], 

Then, observing the overall factor ex^v in eq. (|B2p . we 
can readily drop this term. 
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The 3rd class consists of 10 terms having only one di- 
rectly connected pair. 4 terms being canceled by one 
another, we have the following 6 terms remained; 

+2 TT[g-\dkJd) {du,U-'){d^U)\ . (B3) 

Notice that, when either the band index a or its mo- 
mentum k denotes a filled Bloch state, the diagonalized 
time-ordered Green functions has a pole at the upper half 
plane in the complex w-plane, say, oj — .fc -f iO-|-; 



(B4) 



with tcy.,k — Ea^ki£a,k) = 0. Notc that, in the right hand 
side, we have assumed the infinite life-time for these 
Bloch states. This assumption corresponds to ignoring 
the anti-Hermitian part of the coUisional self-energy, i.e. 
f. 

Substituting eq. (jB3[) into cq. (|B2[) . we then integrate 
over the frequency, such that the summation/integral 
regions over a/k are restricted within the filled Bloch 
states. Using further eq. (|B4p . we then observe that the 
Hall conductivity is indeed characterized by the dual ver- 
sion of [/(I) electromagnetic fields introduced in the text; 



0"A = 



h (27r)" 

S-x=^^^,.m^U)HdkM.Ji.=.^y 
£x-^{[id.U)Hdk.U)-c.c.]^J^^^^^^^, 

with = {dk^f-a,k)- 

APPENDIX C: AMPERE'S LAW 

In a non-interacting Fermi system, the Gauss law solely 
determines the distribution of the U{1) magnetic field. 
When a doubly degeneracy point is formed by two neigh- 
boring energy dispersions, say a-th and [a -)- l)-th band 
at fc = fc™'", 

i^d,a(fc""")=id,a+l(fc'"'"), 

this degeneracy point in a 3-dimensional /c-space becomes 
a source of the spatial divergence of the U{1) magnetic 
field associated with these two bands; 



Vfc • i3" = sign{det V"'}5'^{k - k" 



Pc{k). (CI) 



The sign of the magnetic charge is given by the asymp- 
totic form of the effective 2 by 2 Hamiltonian around this 
degeneracy point; 



(C2) 



doubly degeneracy 1. 

(a) L 




' 3d k-spage 



FIG. 6: (a) Ld,a and Ld,a+i are degenerated at a line in 
the dual (3 + l)-dimensional space, which we dub as dou- 
bly degeneracy line (red curved line). The plane in this 
figure corresponds to a 3-dimensional space, while its com- 
plementary axis is specified by the four-th direction, (b) 
= fc™™(ti;o) plane and u = plane share a "magnetic 
charge" at {ijJ,k) — (tJo, fc'"™('^o))- The relative sign between 
the charge viewed from the former and that from the latter 
is determined by the sign of the x-component of magnetic 
charge current (compare eq. (|C1|) and eq. ([C5}). 



where o-^ stands for the Pauli matrices^. 

As was discussed in this paper, eigenvalues of the La- 
grangian L in Fermi liquid acquire another argument uj 
in addition to the crystal momentum k. The doubly de- 
generacy point in a 3-dimensional space, however, is by 
construction stable against any four-th axis. Thus, along 
this cj-direction, this degeneracy point forms a degener- 
acy line (see Fig. [BJa)); 

LdAk'^'^i^),^) = id,a+i(fc™"(a;),c^), 

Regarding w as a sort of "time", one can therefore say 
that the U{1) magnetic charge in 3-dimensional dual 
space is a conserved quantity. Corresponding to this con- 
servation, we might as well introduce the C/(l) magnetic 
monopole current; 

dk™"^ 

j-;(fc,o.) = sign{ det V''}^-6%k - k'^^^cu)), (C3) 

so that its spatial divergence is balanced by the temporal 
derivative of the magnetic monopole density; 



E 



di. 7™™ 



duPa{k,Uj) 



0. 



Then, in analogy to the Maxwell equation in a real space, 
we can introduce the dual version of "Ampere law" for 
the [/(I) electric field. Namely, the temporal derivative 
of the dual magnetic field and spatial rotation of the dual 
electric field is originated from this magnetic monopole 
current; 



-Jo 



(C4) 



^M—x.y,z 



To see this law in a primitive way, one has only to 
consider the Gauss law in another 3-dimensional space, 
such as kx = constant plane, which crosses the doubly 
degeneracy line mentioned above (see Fig. 6(b)). Namely, 
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our Lu-ky-kz space, say k^ — fc™'°(ti'o) plane, contain a 
doubly degeneracy point at 

(^, fc„ fc,) = (a;o, A:™(^o), fcr'(t^o)). 
Therefore, generically, the doubly degeneracy at /c = 



fc"""(ajo) also becomes a source of the "magnetic field" 
defined in this 3-dimensional space, which is now 
{B", —Ez,Sy)- To determine the sign of the correspond- 
ing "magnetic charge" viewed from this to-ky-kz space, 
expand eq. (|C2p with respect to uj around uj — bj^^, 



^ [co-uo ky-k^'^^LOo) fc. -fcf-(^o) 



where we imposed kx = k^™{LJo) in the right hand side, 
only to discuss the Gauss law in k^ — constant plane. 
Namely, comparing this with cqs. (jCllC2|) . we obtain the 
Gauss law in this LU-ky-kz space; 

d^B^ + dk, + dkj^ = -sign{ det } 

x6{cj-cjo)S{ky - fc™"(^o))'5(fc. - fc,"™K)).(C5) 

By solving luq in favor of k^ , we can actually see that this 
equation is nothing but the x-component of eq. (jC4p . 



APPENDIX D: EXAMPLES 
- ELECTRON-PHONON COUPLING SYSTEMS - 

In this appendix, based on a specific model calcula- 
tions, we will present the crude estimate of the dual ver- 
sion of U{1) electric field. Basically, the electric com- 
ponent arises, in a perturbative calculation, from the 
frequency dependence of the self-energy. Moreover, we 
are interested in the hermitian part of the self-energy. 
This is analytically more difficult to extract than well- 
studied anti-hermitian part. Namely, the latter involves 
an energy-conserving delta function, while the former 
takes over the corresponding principal integrals. As a 
result, a consideration for the case of the Coulomb in- 
teraction becomes quite complex, the frequency depen- 
dence arising only at two loops. We therefore consider 
instead the self-energy due to an electron-phonon inter- 
action, for which we already have a non-trivial frequency 
dependence at one loop. Though we did not discuss 
the electron-phonon interaction explicitly in our deriva- 
tions, the considerations there still apply to this problem. 
Namely, our formulations are clearly free from the specific 
origin of the self-energy, provided that its anti-hermitian 
part can be negligible at low frequency region. 

As a simplest model, we considered the 2D Rashba 
model defined on a square lattice, subjected to an exter- 



^^mm Qf.^ra ^/..mm " 






duj duj duj 

1 


• [^"(^o)] • 


ay 


1 _ 












1 

nal magnetic field along 


the 2;-direction; 





{jm),a 



j ii,v=x,y 

We employ the lattice-regularized model, since the hermi- 
tian part of the self-energy usually depends on the ultra- 
violet cut-off. Namely, a obtained from the Matsubara 
Green function is composed of a principal integral with 
respect to momentum, which shows a ultraviolet logarith- 
mic divergence in a 2D model (see for example eq. (|Dip ). 
The above Hamiltonian, however, can be referred as a 
"Rashba" model, in a sense that, when Fourier trans- 
formed, it reproduces the same energy dispersion as the 
standard (i.e. continuum) Rashba model, around the T 
point; 

Wo = EE^"(^) [^0(fc)]c.,/3 #(fc), 
a,p k 

Ho{k) = —t{coskx -\- cos ky — 2) 1 

+hzTz -h a sin kxTy — a sin kyTx , 

with ^t(fc)^_^^^.e*^- V4(j). 

A finite magnetic field hz lifts the band degeneracy at 
k = (0, 0), (tt, 0), (0, tt) and (tt, tt); 

e± k — —t(coskx 4- cosfc„ — 2) ± |Afe|. 



where |Afc| = y/i^ + a^ (sin fc^ + sin ky). Concomi- 
tantly, the magnetic Bloch wavefunctions for these two 
energy bands also acquire finite z-components of spin in 
the presence of hz] 



.in(^), 





■ -sin(^je~*'^^- " 








. cos (4) _ 
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where 9^ and ipk are defined as follows; 



(cos0i;,sin(?fc) = — - — -(/i^, OA/sin'^ fcj, + sin"^ fc. 



|Afe 



(cos v?fe,sinv3fc) 



sin^ fcr + sin^ fc,. 



:(sin kx, sinfcy). 



We have introduced the magnetic field, because, in a 
time-reversal symmetric system, the a.e.f.-contribution 
to the anomalous Hall effect (AHE) as well as that of the 
dual magnetic field (a.m.f.) always vanishes. Namely, 
the electric field is time-reversally "even" in a T-invariant 
system; 

£°'{k)=£°'{-k), 
while the quasi-particle velocity is T-reversal "odd" ; 

VQ(fc) = -Wa{-k). 

Thus, X VQ,(fc) would vanish after the fc- integration 

in a T-invariant system. 

As for "many-body" effects, we consider an electron- 
phonon coupled Hamiltonian, in which the Einstein 
phonon having constant mass ojq > interacts with this 
Rashba fermion; 

q k,q,a,l3 

The hermiteness clearly requires — g-q in general. We 
take the simplest possible forms for this e-p coupling con- 
stant: cjq ~ 5<j=o, and gq=Q ~ g • 1. Even employing this 
over-simplified coupling form, we already find a substan- 
tial magnitude of the electric field, as shown below. 

The lowest order (imaginary-time) self-energy associ- 
ated with this electron-phonon system is 2nd order in 



s(fc,iw„o = j:^YYH+q){^ 



k+q\ 



q 7=± 



f ni 



nb(t^o) + nf{e^^k+q) _^ "-b(^o) + 1 - "-f (e7,fc-i-q) 

iUrn + t^O — e-y^fc+q ILUm — LOq — ej,k+q 



where nt,{uj) = ^p^^i denotes the Bose distribution func- 
tion. 

At equilibrium, the life-time part of the self-energy 
and hermitian part of the self-energy a are analyti- 
cally continued from this imaginary-time one (see also 
eqs. (ITiTTSi) ): 

Taf3{k,Uj) = iT,af3{k,iuJm = UJ + i\S\) 

-iY.aii{k, iuj,n ^uj ~ i\5\), 

aaf3{k,Uj) = ^{i:af3{k,iuJm ^ LO + i\d\) 

+Y,ap[k,iuj„i =uj -i\5\)]. 



The life-time term clearly becomes zero at T = and 
~ /i. Namely its integral region over the internal line 
is exponentially small: 

(nb(t^o) + «f(e)) • (5(w -I- Wo - e) ~ 0, 
(7ib(c^o) + 1 - ?^f(e)) • 5(t^ - t^o - e) ~ 0, 

when \u} — ^\ <^ ujQ a,t T ^ 0. On the other hand, the 
hermitian part of the self-energy, which is made up of the 
principal integral, remains finite even at T = 0; 



a{k, Lo) 



+ 



P 



I. 



p 



^ dk' 



; }(2^'(°^) 



where we introduced the new integral variable k' = k + q. 

When integrated over k', the off-diagonal elements of 
the 2 by 2 matrix |u^,)(u^, | vanish. Namely, they are 
always odd functions of k'^ or ky, while e^^k' is even; 



sign7 



cos Ok'Tz 



2 2 

-t- sin 9k' cos (fik' Tx - sin 9k' sin ipk' Ty) . 

Accordingly, eq. ljPip becomes diagonal; 



a(fc,o;) = {...}.i + 

^ 7=± 

sign7 cos 9k' dk' 



P- 



a; -I- Wo - ^k',-i (27r)2 
sign7 cos 9k' dk' 



,_^>M w - Wo - efc',^ (27r)2 
= 5o(w).i + 5i(w) f,. (D2) 



The scalar function Si (w) is always a smooth function 
at jw — /z| <C Wq. When this self-energy diagonalized in 
combination with Ho{k), a non-zero derivative of 6*1(0;) 
with respect to w becomes indispensable for a finite a.e.f., 
which we will see below. 

In a 2D model, we generally have 2-component electric 
field and a single-component magnetic field; 

with /i = (kx,ky,Lu). Note that |u^^) diagonalizes the 

L{k,uj) — Ho{k) + <t(w). For any 2-band model, the La- 
grangian can be decomposed in terms of the Pauli ma- 
trices; 



const. 



(D3) 



fi—x,y,z 



where TV now reads (— asinfc^, asinfe^;, hz+Si{Lu)). Then, 
one readily see that the dual fields are identified as the 
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solid angle subtended by this normalized vector N 
A^/|A^|. For example, its spatial component reads 



B7 = 



sign7 



sign7 ■ dNy dN, 
2|7V|3 ' dk^ dky 
sign7 • a^N; 



2\N[ 



cos /Cr COS fc„ 



(D4) 



where the band index 7 = ±. Similarly, the temporal 
components are given as follows; 



sign7 dS\ 



2|iV|3 duj 
sign7 • a 



(dN dN 

■ I X /V, 

dky 



\dkx 
cos fc r sin fc„ 



X N 



2|7V|3 5w 
= sign7 • I I ( cos fca; sin fcj, , — cos fcj, sin k^ ) , 

which clearly indicates that a finite is the essential 
origin of the electric fields. With the lattice constant 
^lattice being explicit, the magnitude of the dual electric 
field estimated on a Fermi surface, i.e. w = is given as 
follows; 



fella 



a-lattii 



7=± 



3 

sign7 



X/, 



3<0 ■ ('^0 

sign7 



] dq 

^>o A,.(<^o + e^,,)2i(2^)2^ 



(D5) 



where Ag denotes the direct band gap at (/-point, mea- 
sured with the transfer integral "t" being an energy unit; 
Aq = 2\Xq\/t. In a same sense, we also made it dimen- 
sionless, the Rashba coupling energy, e-p coupling energy, 
Zeeman energy, phonon energy and band dispersion, like 
a , g , hz , loq and e-y,q = {^-y.q — respectively. In 
terms of these, the other dimensionless function of k, i.e. 
\Nk\, reads; 



\Nk\ 
Si 



a-(sin^ ky + shr fc^) + hl{l + Si)^, 
2 I / sign7 



sign7 



1 dq 
J (271)2- 



,,>o \ ■ i^^o + ej^q) J (27r) 



The expression eq. (jPSp clearly demonstrates that the 
magnitude of the electric field is finite only in the pres- 
ence of the applied magnetic field hz and Rashba coupling 
a. It is also proportional to the dimensionless factor, i.e. 
/, whose value depends on a specific shape of the upper 
band and lower band in /c-space. From its form, however. 



there is no reason that this real-valued factor always has 
to reduce identically to zero. In fact, a simple numerical 

K_y 



COS ky sin k^) , 
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FIG. 7: A distribution of [E^ , £y ) estimated on a Fermi sur- 
face, i.e. uj = jj,. a — 1.0, g = 2.0, hz ~ 0.5, tJo = 0.5 and 
^ = 0. 



estimation shows that / — —0.31 and 5*1 = —0.24, in the 
case of a = 1.0, g = 2.0, hz — 0.5, /i = with ojq — 0.5. 
Observing these numerical factors, the magnitude of the 
a.e.f. is estimated as follows; 



IS" 



10-1 



^lattice 

t ■ 



(D6) 



The distribution of the a.e.f. has a 4 vortex struc- 
ture within the unit cell [— 7r,7r] x [— tt, tt] (see Fig. [T]). 
These vortices reflect the 4 band-crossing points located 
at {kx,ky,uj) — (0, 0, cJi), (±7r, ±7r, wi), (0, ±7r, cji) and 
(±7r,0,wi) with LOi defined as follows; 

hz + Si{uJi) = {). 

Namely, {—Ey,E^,B^) has sources (27r "charges") from 
the former two points, while it has sinks {—2tt "charges") 
at the latter two. 

The many-body correction to the anomalous velocity 
is given by the outer product between the a.e.f. and the 
quasi-particle velocity, i.e., xvc,)^. Thus, taking |va| 
to be aiattico'i, we find it of the order of 10"^ x aj^^^jji^^. The 
bare contribution of the anomalous velocity on a Fermi 
surface can be directly estimated from eq.(|D4|; ~ 
1.0 X a^^j^jj,^. Observing these two quantities, we can 
then insist that the electric-field contribution becomes 
almost at the same order of this bare contribution and 
not negligible even in this oversimplified model. 
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